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We report a helicity analysis of sub-threshold p° production on 2 H, 3 He and 12 C at low photo- 
production energies. The results are indicative of a large longitudinal p° polarization (I = 1, m = 
0) and are consistent with a strong helicity- flip production mechanism. This signature is used 
to extract in-medium p% invariant mass distributions for all three nuclei in a manner which is 
less model-dependent than previous measurements. The results are compared to kinematic and 
phenomenological models of the p° spectral function. The 2 H and 3 He data distributions support 
the role of iV*(1520) excitation in shaping the in-medium p\ invariant mass distribution, while the 
12 C distributions are consistent with quasi-free p° L production. The data support an in- medium 
modification of the pi invariant mass distribution. 

PACS numbers: 13.60.Le, 25.20.Lj, 14.40.Cs, 12.40 Yx 



I. INTRODUCTION 

Of all particles, the p-meson has received the most at- 
tention with regard to medium modifications. Since the 
p° carries the quantum numbers of the conserved vector 
current, its properties are related to chiral symmetry, and 
can be investigated with a variety of models. Most mod- 
els predict a reduction of the renormalized vector meson 
mass in the nuclear medium. A review of the field, up 
to 1999, is given in Ref. O. One of the first models was 
"Brown and Rho scaling" 2] , in which the in- medium p° 
mass is rescaled according to the relation 

m *p M m *N M m *P K m l M ft 
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This relation is based on chiral symmetry and scale in- 
variance, and predicts that the mass of the p° should drop 
by ~ 15% from its free mass value at standard nuclear 
density. This is supported by Lattice QCD calculations 
Q, which suggest that chiral symmetry will be fully re- 
stored at T c > 150 MeV and/or p c > 5 p nuc - 

In addition to mass rescaling, the in-medium shape of 



the p° may be changed by resonant interactions. One 
well-known model of these interactions is that of Rapp, 
Chanfray, and Wambach 4]. Here, the strong coupling 
of the p° with 7r + 7r _ states in the nuclear medium, and 
p°-baryon scattering, lead to a series of "rhosobar" ex- 
citations. This results in an enhancement of the low in- 
variant mass portion of the p° spectral function. This 
model has been applied to the high temperature and high 
density regime of high energy heavy ion collisions with 
great success. This picture has been refined and extended 
by other groups in recent years, such as the relativistic- 
model calculation by Post, Leupold and Mosel These 
and other models will be discussed in more detail later 
in the paper. 

Experimental evidence for in-medium p° mass mod- 
ification has been widespread, but all interpretations 
suffer from significant model uncertainties. For exam- 
ple, CERN dilepton production data from S+Au and 
S+W collisions at 200 GeV/u yield a significant enhance- 
ment at low m ee |6( , indicative of a density-dependent p° 
mass reduction which is consistent with chiral symme- 
try restoration as well as with /3°-medium rescattering. 
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An IUCF 28 Si(p, p') 2S Si polarization transfer experiment 
found an effective isovector NN interaction strength con- 
sistent with m p w 615 MeV/c 2 0, and a KEK 12 
GeV p + A collision experiment yielded e + e~ spectra 
which indicate a significant enhancement for Cu which 
is not present for 12 C, consistent with w 200 MeV/c 2 
p/uj mass shift 8], but these both suffer ambiguities in 
interpretation due to the use of a proton probe. Fras- 
cati -fA total photo-absorption cross sections on C, Al, 
Cu, Sn, Pb at 0.5-2.6 GeV are smaller than expected, 
indicative of increased shadowing. These results are 
best explained in terms of a reduced in- medium mass, 
m p ^610 — 710 MeV/c 2 , which increases the coherence 
length X p = 2k/m 2 p Q However, due to the inclusive 
nature of the experiment, this can only be considered an 
indirect observation of in-medium p° properties. 

Finally, the TAGX Collaboration has investigated the 
3 He(7, 7r + 7r~) reaction in the sub-threshold region. In 
this case, diffractive p° production is suppressed, as the 
reaction must utilize Fermi momentum to produce a p . 
The resulting p 0, s are produced with low boost with re- 
spect to the nuclear medium, so that they decay within 1 
fm of their production point (i.e. still within even a small 
nucleus). These two properties of the sub-threshold re- 
gion reaction enhance the nuclear-medium effect. The 
effects of irN final state interactions (FSI) are minimized 
by the choice of a small nucleus ( 3 He). An energy- 
dependent mass reduction to* w 640 — 680 MeV/c 2 was 
found for 800 < E 1 < 1120 MeV [H EJ E2, with ev- 
idence for even lower p° mass at lower photon energy 
[T^|. Of these, Ref. used an analysis technique sig- 
nificantly different than Refs. [Iol[Tll |. and yielded nearly 
the same m* value. A limitation of these results is the 
model-dependent separation of the p a channel from other 
processes leading to the tt + tt~ state. Experimental data 
which can be interpreted in a less model-dependent fash- 
ion are highly desirable, and this is the goal of the study 
reported here. 

A technique which may allow a cleaner separation of 
the p° — > 7r + 7r~ contribution from competing processes 
has recently been reported by the TAGX Collaboration 
|14|. There, the p° decay angular distribution is recon- 
structed for sub-threshold 2 H and 12 C(7, ir + ir~) reaction 
data, and it is found that when kinematic cuts which 
enhance the relative population of p° — » 7r + 7r~ decay 
compared to competing processes, such as Air produc- 
tion, are applied to the data, a strong cos 2 6** + distribu- 
tion results. The only mechanism compatible with the 
features of the data is the decay of longitudinally polar- 
ized consistent with a strong helicity-flip mechanism 
of p° production. This result forms the motivation for 
this work. Such helicity-flip amplitudes are interesting, 
because they may be related to the p° production mech- 
anism and so assist in such investigations. Here, the em- 
phasis is placed on the use of the longitudinal polariza- 
tion of the p° L to isolate its contribution to the 7r + 7r~ 
data, with the goal being to extract the p\ in-medium 
line shape and compare it to various phenomenological 



models. A preliminary analysis using this technique was 
presented in Ref. (la ] and this work presents our final 
results and conclusions for 2 H, 3 He and 12 C. 

The deuteron comprises an important element in the 
investigation of medium modifications. Its low Fermi mo- 
mentum and nuclear matter density reduce the proba- 
bility of medium modifications and/or nucleonic effects. 
However, the low binding energy and the better defined 
final states for Monte Carlo (MC) simulations also make 
the deuteron an attractive test case for comparisons with 
data and with the free p° line shape. 12 C is a "bench- 
mark" test of nuclear modifications, because of its com- 
bination of nuclear size and density. Unlike the simpler 
targets, the number of nucleons present in 12 C raise po- 
tential complications due to FSI, including absorption. 
The effects of such FSI will be investigated in detail later 
in this paper. 3 He may be expected to represent the 
"nuclear" aspects of 12 C with the simplicity of 2 H, as far 
as FSI and pion absorption corrections are concerned. 
The low relative p°-N and p°— 3 He momenta of this ex- 
periment largely compensate for the short nuclear radius 
disadvantage, compared to more massive nuclei, by in- 
creasing the probability of p° decay within the nuclear 
volume [HI- 

This paper is organized in eight sections. In Sec. II, the 
experiment setup and analysis procedure are reviewed. 
Sec. Ill discusses the MC simulations and presents rep- 
resentative characteristics of the reaction channels con- 
sidered. Sec. IV compares the data distributions with 
the MC simulation results, with particular emphasis on 
the effects of the kinematic cuts on the observed helicity 
angle distributions. In Sec. V, the helicity analysis is 
used to quantify and subtract the proportion of non-p^ 
events surviving the kinematic cuts, yielding the experi- 
mental in-medium p° L invariant mass distributions. Sec. 
VI compares the experimental distributions to a series of 
model calculations, and Sec. VII compares the helicity 
angle analysis result to the previous TAGX Collabora- 
tion results. Finally, in Sec. VIII the discussion and 
conclusions are presented. 



II. EXPERIMENT AND DATA ANALYSIS 

The analysis presented here is based on data taken in 
three separate running periods [H Q Together, 
they comprise data with 2 H, 3 He and 12 C targets. The 
(7, 7r + 7r~) experiments were carried out using the tagged 
photon beam of the 1.3-GeV Electron Synchrotron (ES) 
at the Institute of Nuclear Study at Tokyo (INS) and the 
TAGX magnetic spectrometer. 



A. Tagged photon beam 

The photon beam is produced utilizing the 1.3 GeV 
Tokyo Electron Synchrotron with a duty factor of ~10% 
(lTlflsjl . Figure n displays the tagged photon energy dis- 
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tributions used in this work. As the data for the differ- 
ent targets were obtained at different times, and under 
slightly different experimental conditions, the shapes of 
the resulting tagged photon energy spectra are signifi- 
cantly different for the 3 He and CD? targets. The 3 He 
target data are from Refs. and 0], which used ES 
nominal energies of 800 and f 220 MeV, respectively. The 
CD2 target data were taken with three ES energies of 
1040, 1180 and 1200 MeV, with overlapping tagged pho- 
ton energy distributions. 
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FIG. 1: Tagged photon energy distributions used in this 
work. Data were divided into three E-y bins of 600-800 MeV, 
800-960 MeV, and 960-1120 MeV. 

The data were divided into three tagged photon energy 
bins of 600-800 MeV, 800-960 MeV, and 960-1120 MeV. 
The choice of the limits of these bins was motivated by 
the photon energy thresholds for p° production via the 
quasi-free and non-quasi-free mechanisms in nuclei. On 
1 H, the threshold for production of p° with mass one a 
lower than the nominal 770 MeV/c 2 (i.e. 705 MeV/c 2 ), is 
971 MeV. The division between the mid- and high-energy 
photon bins is placed near here, and so the highest pho- 
ton energy bin is expected to receive a significant contri- 
bution from quasi-free p° production on all three targets. 
On 12 C, the corresponding photo-production threshold is 
723 MeV, while for the p° centroid it is 797 MeV. This 
motivates the division between the low- and mid-energy 
bins to be 800 MeV. In this case, the middle photon en- 
ergy bin may contain contributions from the low-mass 



components of the p° (more than one a below mean) in 
a quasi-free production mechanism, as well as nominal p° 
production from the nucleus as a whole. The mid-energy 
bin also has the best photon luminosity and event statis- 
tics for all three nuclear targets. The lowest photon en- 
ergy bin is deeply sub-threshold and so will be the most 
sensitive to low-mass components of the p° in a nuclear 
environment. However, at these low energies the 7r + 7r~ 
production cross-section is small, and any interpretation 
of the data will be complicated by significant phase-space 
restrictions and non-p contributions. 



B. Targets 

The 3 He data were obtained with the use of a cryogenic 
liguid target of 5 cm diameter and 0.0786 g/cm 3 density 
|19| . This necessitated an empty-target background sub- 
traction procedure, explained in Ref. [ll|. The 2 H and 
12 C data were obtained together, via the use of a solid 
deuterated polyethylene and research grade graphite tar- 
get assembly. The graphite target (0.069 g/cm 2 ) was 
positioned at the nominal center of the TAGX mag- 
netic field, with two deuterated polyethylene targets (to- 
tal thickness 0.618 g/cm 2 ) flanking it upstream (-16.8 
mm) and downstream (+15.8 mm) relative to the pho- 
ton beam. The 2 H data were obtained by subtraction of 
the graphite target data from the CD 2 target data, via 
the procedure discussed in section Hi El 



C. TAGX spectrometer 

The TAGX spectrometer has an acceptance of ir sr 
for charged particles. It consists of a dipole magnet (~5 
kG), drift chambers for tracking, and plastic scintillation 
counters for time-of-flight (TOF) and trigger informa- 
tion. Fig. displays the layout of the TAGX system for 
the 2 H and 12 C measurements. A brief description of the 
detectors follows. An extended description of the TAGX 
spectrometer is found in Ref. [Isj . 

Directly surrounding the target is the inner hodoscope 
(IH), made of two sets of six scintillator counters, one 
on each side of the beam. The IH is used in the trigger, 
as well as in measuring the time of flight (TOF) of the 
outgoing particles [l^, |2£| ■ 

Next are the two straw tube drift chambers (SDC), 
located on opposite sides of the beam. Their operation 
greatly improves the vertex resolution capability of the 
system, which aids in the suppression of experimental 
background. The SDC were designed to preserve the 
TAGX 7r-sr acceptance prior to its installation, to not 
impose extensive modifications of the spectrometer, and 
to not induce significant energy losses to traversing par- 
ticles by keeping its thickness to minimum. The SDC is 
explained in detail in Ref. |2Cj . 

Surrounding the SDC are two semi-circular cylindrical 
drift chambers (CDC), subtending angles from 15° to 
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FIG. 2: Top view of the TAGX spectrometer elements in 
place during the 2 H and 12 C data taking. 



165° on both sides of the beam in the horizontal plane, 
and ±18.3° in the vertical plane. Together with the SDC, 
they are used to determine the planar momentum and in- 
plane emission angle of the traversing charged particles, 
and the vertex position of the trajectory crossings. 

The outer hodoscope (OH) is a set of 33 scintillator ele- 
ments placed after the CDC. Each scintillator is oriented 
vertically, with PMT's attached at top and bottom to 
determine the track angle relative to the median plane. 
The two sets of hodoscopes, IH and OH, measure the 
TOF of the tracked particles. 

Finally, for the 2 H and 12 C data taking a pressurized 
CO2 Cerenkov detector was added to further improve the 
suppression of beam-related electromagnetic (EM) back- 
ground. The gas pressure was slightly higher than atmo- 
sphere, resulting in an electron threshold of 17 MeV/c. 
An aerogel Cerenkov counter with n = 1.030 was located 
at backward angle on beam right. The pion threshold for 
this detector was 570 MeV/c. Because of their proximity 
to the TAGX magnetic field, neither detector was used 
in the on-line trigger decision. 

The remaining components of the TAGX spectrometer 
are four 155 mm x 50 mm x 5 mm scintillator counters, 
with a primary function to veto e + e~ background. These 
veto counters are positioned along the OH arms in the 
median plane, and eliminate charged-particle tracks reg- 
istering within Az = ± 2.5 mm, mostly affecting forward- 
focused e + e~ pairs produced copiously downstream of 
the target, but having a small effect on ir + n~ events. 



D. Event reconstruction 

As our interest lay with 7r + 7r~ production from the 
decay of the p° meson, the experiment trigger was set up 
to record two-charged particle coincidences on opposite 
sides of the beam axis. 

The reconstruction of the TAGX data is briefly re- 
viewed here; a detailed account is in Ref. |18|. In the 
analysis, a lab system of co-ordinates was used, in which 
the x-axis is taken to be along the direction of the photon 
beam and the z-axis is taken along the TAGX magnetic 
field. The trajectory of a charged particle reconstructed 
by the SDC+CDC drift chamber system defined the hor- 
izontal component of the momentum, P xy , the horizontal 
tracked trajectory length, l xy , and the tangential direc- 
tion of the track in the horizontal plane, (p. The result- 
ing planar momentum resolution is given by the relation 
a(P xy )/P xy ~ 0.090 [GeV- 1 ]/^ + 0.12 x 10" 3 

The timing difference of the PMT's mounted on both 
ends of each OH scintillator gave the the ^-component 
of the hit position, which measured the trajectory 
length along the z-axis (l z ), and allowed the magni- 
tude of the momentum in three dimensions to be cal- 
culated. The out-of-plane trajectory length resolution 
was approximately c(/ z ) ~ 1.4 cm. The resulting 
three dimensional momentum resolution is dominated by 
a(P xy )/P xy , above. 

The information from each of the two tracks is com- 
bined to provide the vertex position of each event in the 
x — y plane. The SDC+CDC x traceback resolution is 
1.0 mm (a) [2l|. 

E. Particle Identification 

Fig-Elshows th e observed correlation between the three 
dimensional momentum and the TOF between the IH 
and OH. Particle identification boxes were created as 
shown in the figure to select pions and protons for fur- 
ther analysis. This method does not allow e + (e~) to be 
cleanly separated from tt + (tt~) because of the limitation 
set by the timing resolution. Therefore, in the case of 
the CD2 experiment, TDC information from the forward 
gas Cerenkov detector was used to additionally suppress 
the e + (e~) background. The aerogel Cerenkov detector 
suffered from poor gain due to the TAGX magnetic field, 
and was not used in the analysis. 

Events were only kept if a tt + and a n~ were detected 
on opposite sides of the beam-line, a requirement consis- 
tent with the experiment trigger. Because of the larger 
probability of one or more protons from the 12 C target 
also intercepting the spectrometer, events in which a pro- 
ton was detected in addition to the left-right-going tt + tt~ 
pair were also accepted for the CD2 experiment, only. 
7r + 7r~p events meeting this condition comprised 6% of 
the CD2 target event sample. "Three track" events with 
a 7T + 7r _ pair detected on the same side of the beam-line 
were excluded from further analysis in all cases. 
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FIG. 3: Data distribution of the product of charge and mo- 
mentum versus time of flight through the spectrometer. The 
solid lines indicate the selected regions from which n + , ir~ 
and proton candidates were selected for further analysis. The 
e + e~ background is restricted primarily to the region of small 
TOF and small momentum. 



F. Background Subtraction and Target Separation 

For the 3 He experiment, the cryogenic target was 
housed in a low-mass target cell with 185 fim thick mylar 
walls, a mylar-aluminum laminate super-insulator, and a 
50 /im thick aluminum radiation shield |19| . This neces- 
sitated the use of dedicated empty-target cell data acqui- 
sition runs, for the purpose of background subtraction. 
The result of the 3 He empty-target cell subtraction is 
shown in Ref. [Tlj . 

Dedicated empty-target subtraction runs were not nec- 
essary for the CD2 experiment, as no target container was 
used. In this case, the x — y vertex position provided the 
means to identify whether the event originated from the 
graphite target, the deuterated polyethylene target, or 
elsewhere. Monte Carlo (MC) simulations incorporating 
the known target densities and positions, and taking into 
account the variation of the traceback resolution upon 
the two-track opening angle and the TAGX tracking res- 
olution, were normalized to the data to provide the tar- 
get subtraction parameters to extract the 2 H data from 
the CD2 target yield. Fig. ^compares the x-coordinate 
distribution of the data to the sum of three MC simula- 
tions, in which the 7r + 7r~ production was from cither 2 H, 
12 C, or the atmosphere between the solid targets. The 
inferred contributions of each source to the events within 
each target box shown in Fig. 0|are listed in Table [I] 

The three simulations provide an excellent description 
of the CD2 target regions, and a poorer but still accept- 
able description of the nat C target region. In all three 
cases, the level of agreement is within the 5% systematic 
uncertainty in the MC simulation of the detector, and 
so it was concluded that it was not necessary to add an 
arbitrary background source to describe the target trace- 
back data. After application of a correction to account 
for the variation in TAGX acceptance with rr-coordinate 
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FIG. 4: (Color online) Distribution of the a;-component of 
the vertex position for the CD2 experiment data, and com- 
pared to a series of MC simulations taking into account all 
experiment parameters. The curves are: experimental data 
(solid line), sum of MC simulations (dashed line), simulated 
2 H contribution (dotted line) and air (shaded region). To aid 
the clarity of the figure, the simulated 12 C contribution is not 
shown. The three indicated regions are the limits of the target 
cuts used in the analysis. 



TABLE I: Estimated contribution of the various event sources 
to each of the target boxes, defined in Fig. 01 for the CD2 
experiment. The difference between the sum of the contribu- 
tions and 100% indicates the level of agreement between the 
MC simulations and the experimental data. 



Box # 


Observed Events 


2 H 


12 C Air (1 atm) 


Sum 


1 


38487 


39% 


59% 1% 


99% 


2 


9484 


5% 


85% 5% 


95% 


3 


49873 


38% 


61% 1% 


100% 



(values 1.06 and 0.87 for Boxes 1 and 3 relative to Box 2), 
the Box 2 data were subtracted from the Box 1 and Box 
3 data to yield the 2 H event sample. As air is predom- 
inantly N and O, which will be nearly indistinguishable 
from 12 C in the physics analysis, the decision was made 
to take all of the events within Box 2 as the 12 C event 
sample. 



III. MONTE CARLO SIMULATIONS 

The MC simulations constitute an integral part of the 
analysis, as they provide the best means to take into ac- 
count the effect of the limited experimental acceptance 
upon the observed data distributions. A detailed com- 
parison of the observed data distributions with those 
predicted by the MC simulations will help determine 
the kinematic regions populated by the respective reac- 
tion channels, and how best to separate their contribu- 
tions. Two key observables in this comparison are the 
lab frame pion-pion opening angle (6^) and the missing 
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mass (m miss ). It is necessary for us to investigate the 
effects of analysis cuts on these two variables and any 
possible inter-dependence between them. Finally, it is 
important to show that the cuts suppress non-p° back- 
ground in an understood fashion and that they do not 
induce distortion in the spectra that can imitate the sig- 
nature of a longitudinally polarized p° L . 



A. 7r + 7r production channels considered 

1. p° production 

We hrst discuss the p° production generator. p° pro- 
duction was assumed to occur via a quasi-free "fNp — > 
p°N — > tt + it~ N mechanism, where Nf is the participat- 
ing struck proton with initial Fermi momentum pp and 
the remainder of the nucleus is a sp ectator. Several dif- 
ferent p° line-shapes [22|, I23L [2J, |25|, |26| were considered 
in the simulations. While they all yield very similar re- 
sults, all p° simulations shown here use the line-shape of 
Benayoun et al. [2j| unless specified otherwise. This line- 
shape is based on the I = 1 partial-wave analysis (PWA) 
of e+e- -> tt+tt- data from = 350 to 1000 MeV/c 2 , 
and so is the most authoritative description of the free p° 
low mass tail [53] . The p° then proceeded to decay into a 
7r + 7r~ pair. For the simulations assuming longitudinally 
polarized (1 = 1, m = 0) p production, the decay dis- 
tribution was weighted by cos 2 ^ in the p a rest frame; 
for unpolarized p° production the decay distribution was 
taken to be isotropic. 



2. a production 

Quasi-free cr° — > 7t + it~ production and its decay is 
included as a representative scalar process with Breit- 
Wigner width and centroid as listed in Ref. [2(| [54] . This 
decay is isotropic in the cr° rest frame and it is a vital test 
to investigate whether any cut-induced I = 1, m = 0-like 
signatures appear in the MC simulations. 



3. Baryon resonances 

For the baryon resonances, we follow the work of 
Refs. 0, fill fl3| . where a variety of other reaction 
channels leading to 7r + 7r~ were identified and simulated. 
Thus, quasi-free reactions leading to the production of 
A++TT-, A-7T+, 7V*(1520)tt-, 7V*(1520) -> A++7T" and 



[53] Ref. L'rl by the same authors is based on the same data, and 

yields a nearly identical line-shape. 
[54] The actual values used are m = 800 MeV/c 2 and V = 800 

MeV/c 2 , which are in the middle of the range given for the 

/o(400 - 1200). 



iV*(1520) — > pp° were considered here. The line-shapes 
of the baryon resonances were assumed to be described 
by Breit-Wigner distributions with centroids and widths 
taken from Ref. [2(j, while the shape of the p° from 
iV*(1520) decay uses the PWA shape from [23|. In ad- 
dition, A ++ A~ production was simulated in a — > 
AA(A — 2) mechanism, where A — 2 was a spectator re- 
coiling in some excited state with the appropriate Fermi 
momentum. 

As only quasi-free mechanisms have been simulated, 
m miss cuts must be placed on both the experimental data 
and the MC simulated data to exclude regions populated 
by more complicated reaction mechanisms involving ad- 
ditional nucleons or those producing three or more pions. 
AA is the lowest mass mechanism involving more than 
one target nucleon which leads to the ir + ir~ final state, 
and so it serves as an important test case. Any kine- 
matic cuts which eliminate the AA process should be 
even more effective at eliminating yet more complicated 
multi-nucleon mechanisms. 

The single-nucleon Fermi momentum distributions in 
the quasi-free simulations were chosen appropriate to 
each respective target nucleus. For 2 H, a parametric fit 
to the 2 H(e, e'p) data of Bernheim et al. 27] by A. Saha 
|28j was used. For 3 He, a parameterization from a vari- 
ational calculation using the Argonne potential |29j for 
the 2-body breakup 3 He(e, e'pp)n was used. For 12 C, the 
PWIA parameterizations of J.W. van Orden [3(| for the 
1.P3/2 and lsi/2 shells of 16 was used, and the recoil 
nucleus was placed in the corresponding excited state. 
The remaining nucleon or nucleus, as the case may be, 
was assumed to be a spectator recoiling with equal and 
opposite momentum to the reaction product. 



4- nN final state interactions (FSI) 

For all processes leading to 7r + 7r~ production in the 
600 < Ey < 1120 MeV energy range, the emitted pions 
are in the resonant energy region and pion re-scattering, 
as well as absorption, have large cross sections. Although 
FSI in the form of absorption may play a role in the ex- 
traction of the final p° invariant mass distribution and 
will be discussed in more detail later, re-scattering is of 
particular concern because helicity signatures are used to 
extract the p° L content, to assign confidence levels and to 
define the actual background shapes and strengths to be 
subtracted from the data. While one can argue that for 
2 H and 3 He re-scattering FSI should be a small contri- 
bution, their effect can not be ignored in the case of 12 C. 
For the sub-threshold region considered here, our simula- 
tions indicate a mean ttN CM momentum of 230 MeV/c, 
which yields an approximate 50% probablity that one of 
two 71"^ originating in the center of a uniform density 
nucleus with the radius of 12 C will undergo re-scattering. 

Therefore, a series of MC simulations were made in 
which a pion was randomly selected to scatter with a 
proton, 7r + Nf — * n' + N' , where the pion initial energy 
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and angle are from the production and decay MC gener- 
ator, and the nucleon has appropriate Fermi momentum 
for 12 C. The representative processes for which FSI were 
studied were quasi-free p° production, p° production via 
iV*(1530) -> p°p decay, A++7T - and A++A~. The ttN 
amplitudes are taken from the phase shift analysis pa- 
rameterization of |3lj |. 

The effect of absorption FSI on the p° line-shape was 
also considered. In this case, a parameterization of the 
total pion absorption cross-sections of Ref. [^3 and 
was used to calculate the absorption probability as the 
pion traveled from some random point inside a uniform 
density sphere with radius of 12 C to the exterior. Inelas- 
tic processes such as pion induced pion production were 
not considered since the m rn i SS cut will eliminate them. 

In all cases, the simulated particles were transported 
through the detection apparatus using GEANT 341. The 
TAGX detector performance evaluated in the event re- 
construction process was utilized to accurately simulate 
the TAGX detection efficiency and to reproduce the ex- 
perimental resolution. All experimental thresholds and 
acceptance cuts were applied, and the simulated data 
were analyzed in the same manner as the experimental 
data. 



B. Representative characteristics of the reaction 
channels 

Figure [S] shows representative kinematic distributions 
from four MC simulations with the 3 He target and for 
the 800-960 MeV tagged photon energy range. The ob- 
servables plotted are the missing mass m m i SS , the di-pion 
invariant mass m^, the 7r + 7r~ lab frame opening angle 
6mr, and, of particular importance to this work, the he- 
licity angle #* + , calculated in the di-pion rest frame in 
the following manner: 

1. the longitudinal component p n + of the 7r + momen- 
tum with respect to in the lab frame is formed. 

2. both p + and p + are transformed to the ir + ir~ cen- 
ter of mass frame. 

3. 8*+ is then calculated from 



P+ 

The 3 He nucleus was chosen for a detailed discussion 
on the effects of the cuts on the MC simulations for a 
number of reasons. First, 3 He exhibits features common 
to all three nuclei and most of the conclusions drawn from 
its analysis apply to the other two. Second, there are pre- 
viously published works on the same target, but based on 
a different analysis methodology, that this work can be 
directly compared to. Third, the data quality and statis- 
tics (for the MC simulations to be compared to) are good, 



and this improves the reliability of any conclusions made. 
The Air and AA simulations are selected as representa- 
tive baryonic non-p processes and the a simulation is 
a representative mesonic simulation which will mimic a 
p° in all respects except for the 1 = 1, m = helicity 
signature. 
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FIG. 5: Simulated event distributions for four reactions on 
He, as indicated. Distributions are as-detected in TAGX, i.e. 
with full detection thresholds and experimental acceptance 
conditions. All distributions shown are for the 800 < E 7 < 
960 MeV photon energy bin, but no additional cuts have been 
applied. The four simulations shown were each run for 10 6 
simulated TAGX triggers over the tagged photon range in Fig. 
□ 

The A-7T, cr° and p° simulations make use of a quasi- 
free mechanism, as evident by their very similar missing 
mass distributions, indicating the escape of a single, en- 
ergetic nucleon. The AA simulation involves two ener- 
getic nucleons, and so its distribution peaks at slightly 
higher missing mass. For the remainder of the discus- 
sion, only events with missing mass between 2629 and 
2919 MeV/c 2 were accepted for further analysis. The 
upper limit of range is placed at 110 MeV excitation en- 
ergy of the residual nuclear system and so will be referred 
to as the MM110 missing mass cut. It will eliminate all 
37T production processes. The value of the lower limit 
is not critical, and is placed well below the 3 He ground 
state. 

The quasi-free a and p° processes have similar lab- 
frame opening angle distributions, as each rely on back- 
to-back 7r + 7r _ decay in the meson rest frame. The A-rr 
and AA processes have significantly more strength in the 
W7r < 100° region, and this is characteristic of baryonic 
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decay into non-p channels. The cutoff in the distributions 
for 9^ < 50° is due to the experimental requirement 
of two detected tracks on opposite sides of the photon 
beam-line. 

The bottom row of Fig. [5]shows the respective helicity 
angle distributions of the four processes. Even though the 
(7° — > 7r + 7r~ angular distribution is intrinsically isotropic 
in 0086**+ due to its J = nature, once all experimental 
considerations are taken into account the resulting distri- 
butions are not isotropic. The decay of a longitudinally 
polarized p° is the only process which produces a cos 2 #- 
like signature. The 30% forward-backward asymmetry in 
the simulated p\ distribution, and the fact that the dis- 
tribution maximum is not quite at ±1, are understood 
effects of the TAGX acceptance. 



1. Investigation of the effects of the m m i Sa and 9^^ cuts 

The purpose of this section is to observe the effect the 
missing mass and opening angle cuts have upon the sim- 
ulated cos0* + distributions. The latter form the basis of 
the extraction of the p L yields and distribution from the 
data. 

Fig. shows the same MC simulations as in Fig. 
after application of the MM110 missing mass cut. As 
expected, the removal of the high missing mass events 
has depopulated the low m n + n - region, and so each of 
the respective invariant mass distributions have shifted 
sightly upward in mass. The opening angle distribution 
of the Air simulation has shifted to lower values of 9„„. 
This is a desirable effect, as it will cause the addition of 
the opening angle cut to be more effective in reducing 
the contribution of this very important process from the 
data. The effect of the MM110 cut upon the other open- 
ing angle distributions is much smaller. Finally, we see 
that the effect of this cut upon the helicity angle distri- 
butions is small. 

Helicity angle distributions for the same four simula- 
tions are shown in Fig. {7\ this time with the application 
of a series of opening angle cuts of increasing tightness. 
The O^n > 70° cut is shown for comparison to the previ- 
ously reported results of Ref. ^1 ■ 

One striking observation, by comparing Figs. |S]and[7| 
is the effectiveness of the two cuts in suppressing non- 
p° processes such as Air and A A. The difference in the 
y-axis scale of each figure gives a rough measure of the 
effectiveness of the corresponding cuts. The combination 
of the MM110 missing mass cut and the 8^ > 120° cut 
reduced the Air, AA and a processes by factors of 6.2, 
5.2 and 5.0, respectively, while reducing p\ yield by only 
a factor of 2.4. These factors are in addition to the in- 
herent selectivity of the TAGX spectrometer to coplanar 
ir + ir~ processes such as p° decay. The TAGX limited 
out-of-plane acceptance discriminates against two-step 
ir + ir~ production processes and uncorrelated ir + ir~ pro- 
duction when forming the experiment trigger. This com- 
bination of favorable spectrometer acceptance and anal- 




FIG. 6: Simulated event distributions for four reactions on 
3 He, as indicated. Distributions are as-detected in TAGX, i.e. 
with full detection thresholds and experimental acceptance 
conditions. All distributions shown are with 800 < 
960 MeV and 2629 < m r •>,,,,, m..x- ■ - 

applied. 



£ 7 < 

< 2919 MeV/c 2 (MM110) cuts 




FIG. 7: Simulated event distributions for four reactions on 
3 He, as indicated. Distributions are as in Fig HJ but with the 
addition of an applied opening angle cut, as indicated on the 
right side of the figure. The four simulations shown were each 
run for 10 6 simulated TAGX triggers over the tagged photon 
range in Fig. so the y-axis values indicate the survival 
probability after the various cuts. 
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ysis cuts accounts for the experiment's ability to extract 
a p° signal with small cross-section out of a compara- 
tively large background. Depending on the production 
processes compared, the relative enhancement of p° sig- 
nal can be as much as a factor of ten. 

The next observation to be made from Fig. [7| is that 
only the p® L simulation produces a strong p-wave signa- 
ture. It is clear that an opening angle cut as tight as 
6^ > 120° cannot artificially induce a cos 2 ^ distribu- 
tion in any production process. In the 9 n7r 120° row, 
while the simulated A ++ A _ distribution is featureless, 
one may argue that some biasing of the distributions is 
evident in the A ++ 7r~ and cr° channels. While these two 
distributions do not resemble a p-wave, they do indicate 
that some population at the extremes of the cos6>* + dis- 
tributions can be accounted for by non-p^ decay and that 
care, as well as some additional criteria, may be neces- 
sary to separate p° from non-p° processes. Nonetheless, 
for 9^ > 120°, no background process can create the 
"edge-to-center" ratios that p° L decay exhibits. In fact, 
the population at the center of the distribution is a mea- 
sure of the non-p'j, contributions and can be used as an in- 
put to further analysis. We also investigated whether the 
combination of any two simulated processes can mimic 
the observed p\ distributions. We found that no combi- 
nation of non-p reactions leading to 7r + 7r~ production 
can reproduce the cos6>* + distribution, and the response 
to 6^ and m miss cuts, characteristic of p° L production. 
Any admixture of A ++ A~ together with the other two 
background channels will only move the distribution fur- 
ther away from being p^-like. 

In conclusion, we see that the m m i SS and 9^ > 120° 
cuts together provide a good means of enhancing the pro- 
portion of p° events in the data. This combination of cuts 
does not appear to induce a false cos 2 6> distribution to the 
data, and if used in concert with the "edge-to-center" ra- 
tio of the distribution, can provide a unique means to 
identify the p\ yield. While the 6^ > 120° cut appears 
to be effective, we do not believe it is prudent to cut the 
data more severely than this. Beyond this point, the de- 
graded statistics of the surviving data and the restricted 
kinematics caused by the cut become progressively more 
important factors. 



2. Application of cuts to 2 H and 12 C simulated data 

For the 2 H target nucleus, the effects of the two cuts on 
the cosQ* + distributions for the various reaction processes 
lead to conclusions very similar to those presented for 
3 He. However, given the more complex nature of 12 C, 
two specific investigations are worth exploring in some 
detail. One is based on the fact that this nucleus consists 
of six protons and six neutrons and charge symmetric 
background processes may have differing distributions, 
due to triggering requirements. The other is exploring 
the effects of FSI, an effect that also has its basis on the 
larger number of nucleons present compared to the other 



two nuclei. 

Fig. [S] shows simulations for three non-p° channels as 
well as quasi-free p° production, both with and with- 
out the effect of re-scattering FSI. The upper limit of 
the missing mass cut applied, 10970 < m miss < 11305 
MeV/c 2 , corresponds to 130 MeV excitation of the recoil 
system, and so is termed the MM130 cut. This higher 
upper cut value chosen (as opposed to MM110) is to max- 
imize the data statistics, which are more limited for 12 C 
than the other nuclei, while still remaining below 3ir pro- 
duction. 



1 0970<m miss <1 1 305 (MeV/c 2 ) m mlss and 8„>1 20° 

No FSI w/FSI No FSI w/FSI 
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FIG. 8: (Color online) Simulated event distributions for four 
quasi-free reactions on 12 C, both with and without the effect 
of re-scattering FSI, as indicated. In the case of the with-FSI 
simulations, one of the two pions was randomly selected to 
undergo np scattering, weighted according to the phase shift 
analysis parameterization of Ref. |3l|. All distributions are 
as-detected in TAGX (i.e. detection thresholds and exper- 
imental acceptance included) and are for 800 < E 1 < 960 
MeV and 10970 < m miss < 11305 MeV/c 2 (MM130). In ad- 
dition, the opening angle cut 9 W7r > 120° has been applied to 
the two rightmost columns. The four simulations shown were 
each run for 10 6 simulated TAGX triggers over the tagged 
photon range in Fig. so the y-axis values indicate the sur- 
vival probability after the various cuts. 

With only the MM130 cut applied, all non-p° pro- 
cesses exhibit distributions which peak near the center 
of the cos6?* + range. With six protons and six neu- 
trons present, the charge symmetric A ++ 7r~ and A _ 7r + 
reactions should contribute at approximately the same 
rates. However, the inclusion of re-scattering FSI and 
the slightly different TAGX acceptances for tt + and ir~ 
may induce a difference between the charge symmetric 
distributions. Nonetheless, no artificial p-wave behavior 
is seen as a result of the MM130 cut. This cut is also 
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very effective in suppressing AA processes. 

For the non-p processes, the application of the 9^ > 
120° cut results in the suppression of yield in the re- 
gions between -0.7 and 0.7, and the remaining popula- 
tion about cos6>* + = is a measure of the total surviv- 

77 __ 
ing background. All of the simulations in Fig. [H] were 

run for the same number of simulated TAGX triggers, so 
the poor statistics for the FSI simulations indicate their 
low survival probability after the cuts. The mean cut 
survival probability of FSI processes is less than 25% of 
that of non-FSI processes, enhancing the expected rela- 
tive proportion of non-FSI to FSI events in the 12 C data 
by a factor of four. It is interesting to note that the Att 
channels exhibit distributions that may be fitted by a 
function, such as A + B cos9* , + C cos 2 6** + , that has 
a p-wave component. It is also worth noting that the 
"skewness" observed for the A ++ 7r~ reaction (forward- 
backward asymmetry) will give different values for the 
three coefficients than that expected for the p° L decay 
process. These aspects of the analysis will be pursued 
for all three nuclei in a later section of this work. 
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FIG. 9: 3 He data distributions for 800 
In addition, the indicated opening angle 
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IV. DATA DISTRIBUTIONS AND 
COMPARISONS WITH MC SIMULATIONS 

In this section, the effects of the two cuts, m m i SS and 
Q tttt , on the data distributions are investigated and de- 
tailed comparisons with the MC simulations are made. 
The objective is to establish the consistency of the con- 
clusions derived from the application of the cuts on the 
data and the simulations. In subsequent sections in this 
paper, this will lead to the identification of a p° L com- 
ponent in the experimental data and its extraction and 
separation from the non-p° background. 

A. The effects of the cuts upon the 3 He data 

The effect of the 6 V7r cut upon the 3 He data distribu- 
tions for m m i ss and other variables is shown in Fig. [5] 
Four successive 9 nn cut values of 70°, 90°, 110° and 120° 
have been applied. The effect of the cut upon m miss is 
small. There is a shift in the mean of the distribution 
of approximately 30 MeV/c 2 as the cut is changed from 
70° to 120°. The low missing mass region is most af- 
fected, due to the elimination of two pion events with 
low relative momenta. The effect on the invariant mass, 
TOtttt, is also small. Neither the peak nor the shape of 
the distributions change appreciably for the four values 
of 9 n7r cut. The low mass tail of the distribution, below 
500 MeV/c 2 , is affected a little more greatly, resulting 
in a slightly more symmetric shape for the distribution. 
Generally, invariant mass distributions for such multi- 
channel (inclusive) processes are not sensitive variables 
on their own. This is because they represent phase-space 
and acceptance limitations compounded with broad and 
overlapping resonances. 



The effect of the MM110 missing mass cut, in addition 
to the opening angle cut, is shown in Fig. EH A com- 
parison of Figs. l9l and ITUI reveals that the application of 
the m m i ss cut has a larger effect on m TW than the 6 V1T 
cut had alone. By eliminating three pion production and 
suppressing the emission of energetic nucleons, the miss- 
ing mass cut helps to isolate two energetic pion events 
adding up to larger invariant mass. As an example, the 
peak of to^tt corresponding to the 120° cut has shifted 
from approximately 580 MeV/c 2 in Fig. |^lto approxi- 
mately 665 MeV/c 2 in Fig. Similar behavior was 
observed in 0] for 2 H. These shifts become a significant 
factor in the discussion on the p\ invariant mass distri- 
bution and comparisons with theoretical models in the 
relevant sections of this work. 
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FIG. 10: 3 He data distributions for 800 < E 1 < 960 MeV and 
2629 < m miss < 2919 MeV/c 2 (MM110 cut). In addition, a 
series of opening angle 6^^ cuts have been applied to the data, 
as indicated. 
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The effect of the opening angle cut upon cos0* + is quite 
dramatic. The effect of the additional MM110 cut has 
been to make the resulting cos 2 # distribution more sym- 
metric. It is instructive to compare the corresponding 
cut plots in Figs. and ^] The MC simulations show 
cos6>* + distributions for all background processes that are 
radically different from the data. Furthermore, no com- 
bination of Air, AA and a processes without p° L content 
can reproduce the shape and change of population of the 
data in going from the 70° to the 120° cut. The "edge-to- 
center" ratio is an additional piece of information, as the 
MC simulations indicate that the bulk of the background 
processes are located in the "center" region of the cos#* + 
distribution. Thus, this ratio is a measure of the p\ to 
non-p° contributions to the observed data yield, and a 
comparison of Figs. |§1 and ^| indicates that the effect 
of the m m i ss cut has been to remove additional non-p^, 
processes, such as 3ir production. 

In conclusion, the two kinematic cuts, applied to both 
data and MC simulations, have been shown to be largely 
uncorrelated, each selecting different features of the data- 
set. They help to establish that the only self-consistent 
process identified as the source of the p-wave-like dis- 
tribution in cos#* + is the decay of a longitudinally po- 
larized p\. By comparing the 'edge-to-center' ratios of 
the distributions in Figs. [7] and ^3 it is clear that the 
non-p'j, background under the 120° cut is minimal. This 
is established by the MC simulations for the Air, AA 
and cr° processes and the small number of data events at 
cos6>* + = 0. Since this "pedestal" is an accumulation of 
surviving events from all background processes, it estab- 
lishes an absolute measure of the contribution of these 
processes to the data. 

One question that arises as a result of the comparison 
of the simulations and the data is that of helicity conserv- 
ing pj, or unpolarized p° content in the data and their 
signatures in the data and simulations. Both unpolarized 
p° and transversely polarized p? events, the latter with a 
distribution peaking at zero values of cos^, would add to 
the "pedestal" at the center of the distribution. As such, 
it is clear that there is little or no evidence for such po- 
larization states of p° in the data sample. The argument 
against p^ signatures is particularly strong because such 
a helicity state would otherwise share all other responses 
to 9 V7T and m miss cuts with p\. In this case, increasing 
the cut limit would result in a central enhancement, 
an expectation completely in variance with the data. 



B. The effects of the cuts upon the 2 H data 



The 2 H data distributions respond to the 9^ cuts in a 
manner similar to those in Fig. except that the m m i 8S 
distributions show a peak at the deuteron mass with a 
narrower distribution than the 3 He distributions, and the 
m^^ distributions have slightly higher centroids, a man- 
ifestation of the smaller phase-space available and the 
more well-defined final state. Application of the missing 



mass MM130 cut results in the data distributions of Fig. 
1111 The effect of the m miss cut on is similar to that 
already shown on 3 He. The low invariant mass regions 
below 500 MeV/c 2 are essentially eliminated; the means 
of the distributions are approximately 25 MeV/c 2 below 
their corresponding distributions in Fig. 1101 This is con- 
sistent with the lower center-of-mass energies available 
for pion production due to the lower Fermi momenta of 
the struck nucleons. 
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FIG. 11: 2 H data distributions for 800 < E-y < 960 MeV and 
1745 < m miss < 2006 MeV/c 2 (MM130 cut). In addition, 
a series of opening angle nn cuts have been applied to the 
data, as indicated. 



The MM130 cut does not have the same quality of 
improvement in the p-wave-like behavior of the cos#* + 
data as it did for 3 He. This is partly due to the fact 
that the deuterium data are the result of a subtraction, 
CD 2 - 12 C, and partly due to the MM130 cut itself. This 
cut is perhaps too generous for such a loosely bound nu- 
cleus, and will be tightened in subsequent sections of this 
work. Nevertheless, the "edge-to-center" ratios of cos6** + 
distribution for the 120° cut are good. 



C. The effects of the cuts upon the 12 C data 



Fig. H2lshows the effect of the 6^ cuts on the 12 C data. 
The m m i ss and m„ distributions are similar to the two 
lighter nuclei, but with broader, smoother distributions, 
reflecting the larger phase-space available to the pions 
due to the higher Fermi momentum in this nucleus. The 
overall effect on cos6>* , is also similar for all three nuclei. 

7T > 

One minor difference observed in this figure is that the 
loss of acceptance observed at the extremes of the distri- 
bution is even less pronounced for 12 C compared to 3 He. 
This is consistent with the MC simulations. The p-wave 
signature is clearly observed as a result of the opening 
angle cut. 

Application of the MM130 cut results in the distri- 
butions in Fig. 1131 The reduction in the number of 
surviving events is dramatic, compared to the other two 
nuclei. Application of the same cut on 2 H resulted in a 
survival fraction of 54% for the 8^ > 120° cut. In the 
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FIG. 12: 12 C data distributions for 800 < E-, < 960 MeV. 
The indicated opening angle 6 nn cuts have been applied to 
the data. 



3 He case, the surviving fraction is 34% while in this case 
the fraction is 19%. This severely restricts the surviving 
statistical precision. However, comparing the cos#* + dis- 
tributions for the 90°, 110° and 120° cuts in Figs, feland 
1131 the MM130 cut reduces the non-p-wave background 
in a very effective manner. 
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that the 008*6^+ distributions for the three photon en- 
ergy bins of 600-800, 800-960 and 960-1120 MeV are also 
very similar. A sample plot of the 12 C data for the full 
photon energy range is shown in Fig. 1141 Application 
of the O^ir > 120° cut removes most of the yield near 
cos0* + =O, resulting in the p- wave- like distributions in 
the second column. Application of only the MM130 cut 
eliminates events across the distributions, without alter- 
ing their shapes. The final column displays the final dis- 
tributions after the application of both cuts. The sur- 
viving events are consistent with processes dominated by 
p-wave-like distributions. The two-step processes anal- 
ogous to those simulated in Fig. cannot be the main 
contributors, as they have different responses to the cuts 
than the foreground channels and the data respond to 
the cuts as the simulations predict. There is little evi- 
dence that background processes dominate the data-set 
after application of the two kinematic cuts, each with 
their own unique effects on both the data and the MC 
simulations. In fact, a direct comparison of Figs. |S|and 
^3 indicates that the total background is a small fraction 
of the total event sample. However, the exact contribu- 
tion of the background in the cos#* + distributions needs 
to be determined in a quantitative manner, and this is 
the subject of the next section. 
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FIG. 13: 12 C data distributions for 800 < £ 7 < 960 MeV and 
10970 < m mias < 11305 MeV/c 2 (MM130 cut). In addition, 
a series of opening angle 6 nw cuts have been applied to the 
data, as indicated. 



D. Conclusions from the comparison of the data to 
the simulations 

Overall, the effect of the two cuts on all the variables 
investigated and for all three different nuclei do not result 
in any significant differences in the distributions among 
the three target nuclei. It should also be mentioned 



FIG. 14: Sample C data distributions for the three tagged 
photon energy bins, as indicated. 



V. INVESTIGATION OF NON-p£ 
BACKGROUND VIA HELICITY ANALYSIS 

The purpose of this section is to investigate the relative 
contributions of the various 7r + 7r~ production processes 
to the observed cos#* + distributions, to subtract the 
contributions of the background processes in as model- 
independent manner as possible, and in so-doing to ob- 
tain a better measure of the observed p^ L distribution ver- 
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sus invariant mass for each of the target nuclei. In order 
to gain more insight into the nature of the non-p^ back- 
ground remaining after the missing mass and opening 
angle cuts, both the experiment and the MC-simulated 
cos#* , distributions were fit with a function that con- 
tains a p-wave term, a s-wave term and an interference 
term. In this manner, all of the elements of the distribu- 
tions, namely a uniform event distribution, a symmetric 
p-wave distribution about 8* + = 90° and an asymmetric 
interference term are represented in the fit. The fitted 
function is 

Events = A + Bcos0*+ + Ccos 2 6>;+ . (1) 

This function is not intended to provide a perfect fit to 
the observed distributions, but rather to provide a means 
to quantify the relevant features of the observed distri- 
butions. 

If the data show high fractions of p-wave content, it 
means that the coefficient A, which is the "pedestal" at 
cos#* + = 0, is very small. This can either mean that the 
background is negligible and the data are almost purely 
p\ in content, or that there may be significant back- 
ground, but it is almost purely p-wave-like due to the 
cuts imposed. None of the background processes in Fig. 
[7| appear to reflect such a possibility. 

In order to extract the in-medium p° L invariant mass 
distribution, one needs to investigate the data and MC 
simulations as a function of m wn and possibly extract 
different confidence levels for different m„ values. Both 
the experimental and MC-simulated data were binned 
according to 7r + 7r~ invariant mass, and a cos#* + distri- 
bution formed for each m„ bin. Each invariant mass 
bin was 90 MeV/c 2 wide, staggered in 30 MeV/c 2 incre- 
ments. Thus, the first bin was 360< <450 MeV/c 2 , 
the second was 390< <480 MeV/c 2 , and so on. 

This binning procedure was necessary to ensure suffi- 
cient statistics per bin to allow inferences to be formed, 
while preserving sufficient invariant mass resolution to 
make these inferences interesting and useful. The three 
coefficients A, B, C were determined for each bin for 
both the data and the various MC-simulated processes, 
and tabulated according to the mean value of the 
respective population in each bin. Thus, although the 
data and simulations used the same binning scheme, the 
TO W7r tabulated values for each differ slightly according 
to the distributions of the populations of each. Checks 
were made to ensure that the resulting parameter distri- 
butions versus were not sensitive to the bin size and 
stagger offset used. 

If the non-p° processes had cos6>* + distributions de- 
scribed by only the flat A and skewed B terms of the 
fitting function, and if the p\ process had a pure cos 2 6** + 
distribution with only the C coefficient being non-zero, 
then a plot of C versus invariant mass would give a pre- 
cise measure of the observed p° L distribution for each of 
the target nuclei. To first order, this is in fact a good ap- 
proximation. However, as we have seen in section Al! B II 
the TAGX acceptance modifies the expected p\ distribu- 



tion from a pure cos 2 #* + function, and the background 
processes contribute somewhat to C as well as to A and 
B. 

If one integrates equation^Jfrom cos#* + =-l to +1, one 
obtains 2A+^C. The ratio 3 ^ +c thus provides a relative 
measure of the cos 2 #* + component in the observed distri- 
bution. The skewness coefficient, B, does not contribute 
to the integral as this term is odd in cos6** + . However, 
as we have seen in Fig. |UJ certain background processes 
have significantly skewed distributions, so the value of B 
could assist in the discrimination between p° L and non- 
p\ processes. This conclusion was confirmed after further 
analysis, and so a p-wave shape ratio 

^- wave ~ 3\A\ + I\B\ + \C\ 

was formed for further study. This weighs the skewed 
coefficient as if the absolute value of the skewness term 
had been taken in the integral. The absolute values in the 
denominator ensure that cancellation between opposite 
sign coefficients cannot artificially boost the value of the 
ratio. 



A. Comparison with data 

Fig. displays the observed p-wave shape ratio for 
the 3 He data and several MC simulations. The two p\ 
producing channels considered are quasi-free p\ produc- 
tion and JV*(1520) — > p°N decay, while the two "bench- 
marks" for non-p° di-pion production are the quasi-free 
a and An processes. The elimination of the AA pro- 
cess from further consideration is based on the observa- 
tion that it has no p-wave-like structure independent of 
any cuts applied, and that the proportion of AA events 
surviving the cuts is small (see Figs. 0and|HJ. For the 
quasi-free p° L simulation, R p - wave deviates from 1.0 be- 
cause of the finite TAGX acceptance, primarily due to 
the inability to detect pions emitted between and 15° 
in the lab frame. As expected, the V*(1520) -> p°N 
simulation yields similar results. The ratios for the two 
background processes, however, show much greater vari- 
ation, as a result of the m m i SS and 9^ cuts. For example, 
at 960 < E 1 < 1120 MeV the a and An channels display 
high p-wave ratios between 600 < < 700 MeV/c 2 , 
but lower ratios in most other regions. 

For 800-1120 MeV, the data have p-wave ratios which 
are mid-way between the p^-producing and background 
simulations. This means that after the application of the 
Timiss and cuts, the data still contain some contribu- 
tion from non-p^ processes. The effect of the background 
processes is to dilute the p-wave shape ratio observed 
for the data. For example, between 800 and 960 MeV 
tagged photon energy, a less than 10% contribution of 
the a -0 process to the remaining data would be sufficient 
to account for the degradation of the ratio from that of 
the pure-p'j, process to that observed for the data. At 
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FIG. 15: (Color online) p-wave shape ratios for the 3 He data 
[■], quasi- free p\ simulation [long-dash line], JV*(1520) —* pp° 
simulation [short-dash line], A ++ 7r~ simulation [dash-dot 
line] and a simulation [dotted line] , for the three tagged pho- 
ton energy bins, as indicated. The solid lines are normalized 
fits of the simulations to the R p - wa ve data, as described in the 
text. The regions of largest discrepancy between fit and data 
are generally on the tails of the event distributions. All distri- 
butions shown are with 9 W7T > 120° and 2629 < m m iss < 2919 
MeV/c 2 (MM110) cuts applied. 



600-800 MeV, however, the degradation of the observed 
ratio is more severe, and a more significant contribution 
by background processes to the remaining data after the 
cuts would be necessary to account for the observed ratio. 
For the higher £L region of 960-1120 MeV, the picture 
is even more complicated. Within the data error bars, 
the 600 < m„ < 740 MeV/c 2 region has a high p-wave 
shape ratio, but at the same time, the background chan- 
nels (especially the a channel) also have high ratios. 
Thus, the level of confidence in this region is not as high 
as that of the mid-energy region. 

Thus, the confidence in the accurate extraction of the 
p° L component, and therefore its invariant mass distribu- 
tion, depends on both E-y and m^. In the mid-i? 7 range, 
where the most and highest quality data have been ob- 
tained, the degradation in the observed p-wave ratio ap- 
pears to be small, and so p\ events may be extracted 
with high confidence for values between 550 and 



730 MeV/c 2 . The confidence in the low E 1 region, which 
is deeply sub-threshold, is not as high, while the situa- 
tion at the high end of the tagged photon energy range 
is similar. 

Fig. displays the p-wave shape ratio for the H 
data and the equivalent four MC simulations. In this 
case, the imposed missing mass cut is at 90 MeV excita- 
tion energy equivalent (MM090), in order to improve the 
foreground to background ratio of the data. The open- 
ing angle cut remains at 120°. The 2 H data exhibit a 
p-wave ratio which is diluted in comparison to the pure- 
p\ channels, which indicates that a portion of the data 
remaining after the cuts is due to background processes. 
This is particularly true for the low-energy panel, where 
the 550-620 MeV/c 2 invariant mass data have a ~ 40% 
ratio, while the pure-p^ channels have ratios in excess of 
80%. In this case, the dilution of the ratio could be ex- 
plained by attributing approximately 40% of the data to 
background processes, while the background proportion 
would be much less for the other two energy bins. 
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FIG. 16: (Color online) p-wave shape ratios for the 2 H data 
[■], quasi-free p° L simulation [long-dash line], iV*(1520) — > pp° 
simulation [short-dash line], A ++ ir~ simulation [dash-dot 
line] and a simulation [dotted line] , for the three tagged pho- 
ton energy bins, as indicated. The solid lines are normalized 
fits of the simulations to the Rp- wa ve data, as described in 
the text. All distributions shown are with 8 n7r > 120° and 
1745 < m miss < 1966 MeV/c 2 (MM090) cuts applied. 
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Finally, a measure of the reliability of the helicity anal- 
ysis for 12 C can be pursued with the same methodology 
in Fig. El The data display larger p-wave shape ra- 
tios than either of the 2 H or 3 He data, consistent with a 
smaller proportion of background reactions surviving the 
two cuts. 

Simulations both with and without the effects of re- 
scattering FSI are also displayed in Fig. El 111 the 600- 
800 McV photon energy bin, the data and FSI-modificd 
p° L processes are in excellent agreement for all values of 
TOtttt. This is likely a fortuitous agreement, as it would 
indicate a 0% contribution due to background processes 
in this energy region, an expectation at variance with 
our findings for the two lighter nuclei. Moreover, in the 
800-960 MeV energy range, the data display a p-wave 
shape ratio which is larger than that expected for the 
FSI-modified p\ processes. To be in agreement with 
the data ratios, a negative number of background pro- 
cess events would have to be subtracted in this region, 
which is unreasonable. 

Simulations without FSI are shown in the left panels of 
Fig. 1171 There, the data show large p-wave ratios which 
are consistent with or degraded slightly from those ex- 
pected for pure p° L processes. Our expectation from Fig. 
03 is that a large fraction of the re-scattering FSI will be 
removed by the analysis cuts. If one assumes that 100% 
of FSI processes have been eliminated, the non-p 1 ^ contri- 
bution to the remaining 600-800 MeV and 800-960 MeV 
tagged photon bin yield is still less than 10%. A more 
likely scenario is that some FSI contribution remains af- 
ter the cuts, but since the re-scattering processes do not 
significantly affect the p-wave shape ratio, their contribu- 
tions are relatively difficult to distinguish. In this case, 
however, the assumption of no re-scattering FSI actually 
leads to the subtraction of a greater amount of non-p^ 
background, and so is the more conservative of the two 
choices. The 960-1120 MeV energy bin is less definitive, 
due to poorer statistics. Even here, the data exhibit p- 
wave ratios consistent with p\ dominance, although the 
background contribution is likely higher than in the two 
lower tagged photon energy bins. 

It is interesting to note that among the three nuclei, 
the more massive and complex target exhibits the clean- 
est p-wave signature overall. This analysis simply under- 
scores the observations made by comparing Figs. I10II11I 
and 1131 which are m^-integrated distributions for the 
data surviving the cuts from each nucleus. It appears 
that the missing mass cut is very effective in eliminating 
background, moreso in the case of 12 C. 

The response of the MC simulations to the two cuts 
provides a good description of the response exhibited by 
the experimental data. While all three nuclei seem to 
be dominated by the same production mechanisms, as 
indicated by the distributions of several important ob- 
servables, and all three show clear signatures of p-wave 
attributed to p a L production and decay, the relative back- 
ground proportion is different for the three nuclei. It is 
important to emphasize here that the information ob- 



tained for the background channels are relative ratios, 
since absolute cross sections are not known under these 
conditions. However, they are valuable inputs to the ro- 
bustness of the conclusion that p° L events are correctly 
identified for specific regions of photon energy and in- 
variant mass. They also provide an additional, and inde- 
pendent, tool to compare the results of this analysis with 
previous publications on the subject |l0llllLll3l | by form- 
ing the foundation upon which the extraction of the p° L 
invariant mass distribution is accomplished. For all three 
nuclei, the 800-960 MeV tagged energy regions have the 
best statistics and cleanest p-wave signatures associated 
with 1 = 1, m = quantum states. This is also the en- 
ergy region of increasing, but still sub-threshold, p° pro- 
duction. Thus, this is the region where one expects non- 
coherent production signatures to manifest themselves, if 
they are to be present at all. Finally, the analysis and FSI 
simulations shown here have not altered the conclusions 
in Ref. 0], that no background processes have been 
identified which can account for or explain the helicity 
signatures observed in the data. 



B. Non-pi background subtraction 

Based on the analysis of the p-wave shape ratios of the 
different contributing processes, the portion due to non- 
p L background can be estimated and subtracted. The 
basic procedure was that relative contributions of the 
foreground processes (TV* (1520) — » pp° and quasi-free 
p° L ) and background processes (A ++ 7r~ and a ) were fit 
to the R p - waV e distributions of the data, yielding two 
relative normalization factors r\ p for the foreground and 
(1 — rj p ) for the background, respectively. Our estimate 
of the in-medium p° L invariant mass distribution is then 
obtained from 

p° L (m^,E 7 ) Distribution = C da ta - (1 - r) p )C back ground- 

A single normalization factor was thus obtained for each 
tagged photon energy bin; the values are tabulated in 
Table [D] To minimize model sensitivity, the A, B. C 
coefficients for the two foreground processes, quasi-free 
p° L and TV* (1520) — > pp°, and for the two assumed back- 
ground processes, quasi-free A ++ 7r~ and a , were aver- 
aged together prior to fitting the R p - wave distributions. 

TABLE II: Relative proportion of foreground processes 
(quasi-free pi and iV*(1520) — > ppi) m the data, after miss- 
ing mass and opening angle cuts applied, as determined from 
fitting the R p - ma ve distributions of the data. 



Nucleus 


Tagged photon energy bin (MeV) 




600-800 800-960 960-1120 


2 H 


0.66 ± 0.06 0.75 ± 0.04 0.87 ± 0.06 


3 He 


0.25 ± 0.11 0.91 ± 0.02 0.83 ± 0.05 


12 C 


0.90 ±0.01 0.98 ±0.02 0.96 ±0.14 
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FIG. 17: (Color online) P-wave shape ratios for the 12 C data [■], quasi-free p\ simulation [long-dash line], A r *(1520) — » pp° 
simulation [short-dash line], A ++ 7r~ simulation [dash-dot line] and a simulation [dotted line], for the three tagged photon 
energy bins, as indicated. The panels on the left assume no re-scattering final state interactions, while those on the right assume 
that one of the two detected pions has re-scattered. The solid lines are normalized fits of the simulations to the Rp-wave data, 
as described in the text. All distributions shown are with 0™ > 120° and 10970 < m mias < 11305 MeV/c 2 (MM130) cuts 
applied. 



The solid lines in Figs. 1 1 51 1 1 61 and 1 1 71 indicate the fits 
to the R p - wave distributions for the three target nuclei. 
It is clearly seen that the MC simulation sums do not 
perfectly describe all of the features of the data. How- 
ever, the p-wave shape ratios are relative ratios, indepen- 
dent of the actual yield of the data within each invari- 
ant mass bin, and the discrepancies between the fits of 
the simulations and the data are largely confined to the 
tails of the invariant mass distributions. To avoid unwar- 
ranted sensitivity to these regions, the R p - waV e fits were 
weighted according to the data yield within each invari- 
ant mass bin. In many cases, the fit to the data could 
have been improved by choosing one of the background 
or foreground processes alone in the fit and excluding the 
other, rather than averaging them. If this was done, the 
obtained normalization factors would not vary outside 
the uncertainties listed in Table ITT1 Nonetheless, which 
background or foreground process to choose to obtain the 
best fit would depend on target nucleus and photon en- 



ergy bin, and so the averaging procedure was maintained 
to ensure uniformity of approach for all of the data dis- 
tributions. 

Fig. El shows the cos 2 #* + fit coefficient 'C ver- 
sus invariant mass for all three nuclei. The shaded re- 
gions are for the estimated background contributions to 
the cos 2 #* + coefficient distributions resulting from the 
Rp-wave analysis. With the exception of the 600-800 
MeV 3 He distribution, the estimated background con- 
tribution is nearly negligible for 12 C, a bit larger for 
3 He, and larger yet for 2 H. This systematic behavior is 
despite the tightest missing mass cut being applied to 
the 2 H data, and the most generous cut applied to the 
12 C data. In most cases, the shaded contributions to 
be subtracted are relatively small, quantifying the earlier 
discussion that the data display large p-wave signatures 
which are consistent with p° decay [l4| . 

Regarding the 600-800 MeV distribution for 3 He, one 
should note from Fig. ^ that the photon energy mean 
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FIG. 18: (Color online) cos 2 #* + fit coefficient 'C of the data [U] and the estimated non-p^ background contribution [shaded 
region], based on the fit to the R p - wa ve distribution. 



and distribution for this bin are markedly different than 
for 2 H and 12 C, and so substantial differences are to be 
expected. Because of the 60 MeV lower mean photon en- 
ergy, the 7r + 7r~ cross-section on 3 He is reduced, and the 
data have poor statistical confidence. The large data er- 
ror bars preclude any meaningful conclusions to be drawn 
with this analysis technique. Even with these limita- 
tions, the value of the cos 2 6>* + coefficient for m„ > 575 
MeV/c 2 indicates the presence of some p\ contribution to 
the data. The maximum R p - wa ve values of the 600-800 
MeV 2 H data are nearly double those of the background 
processes, and this can be accommodated by assuming a 
foreground/background proportion of 2:1 fTable ITT}) . The 
background contribution to be subtracted (shaded region 
of the left panel of Fig. I16f) has nearly the same invari- 



ant mass distribution as the data, and so the resulting 
distribution is not particularly sensitive to errors in the 
subtracted proportion normalization. The 600-800 MeV 
12 C data are fitted quite well by the sum of the MC sim- 
ulations, and this is reflected in a substantially smaller 
uncertainty in rj p . 

As has been already discussed, the 800-960 MeV region 
is one of good statistical precision and one of high con- 
fidence in the p-wave shape analysis. This is reflected in 
the fits of the simulations to the data, and the low uncer- 
tainties in the r\ p normalization factors. In all cases, the 
portion to be subtracted is approximately centered be- 
neath the data, and so the background subtraction will 
appreciably alter neither the centroid nor the width of 
the extracted in-medium p\ invariant mass distribution. 
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While the 2 H and 3 He data are fit quite well, the 12 C data 
display a R p - waV e ratio which is slightly larger than any 
of the MC predictions. Recall that for the simulations 
including FSI, this discrepancy was even larger. As the 
maximum R p - wave value of the iV*(1520) — » pp° L simu- 
lation equals that of the data, just at a different invari- 
ant mass value, the discrepancy could be in-part due to 
the free p° line-shape assumed in the simulations. The 
data indicate a very small non-/?° proportion to be sub- 
tracted, and this is reflected in the near-unity value of 
the obtained r\ p factor. 

The quality of the fits to 960-1120 MeV energy bin data 
are poorer than at 800-960 MeV. At this higher energy 
p° — » 7r + 7r~ decay is more likely to result in missing mass 
and opening angle values which do not pass the physics 
analysis cuts, resulting in the poorer statistics of the sur- 
viving data sample. The number of incident tagged pho- 
tons for this bin also is smaller, especially for the CD2 
experiment. With the exception of the extreme tails of 



the distribution, the fit to the H R„ 



distribution 



generally follows the trend of the data, given the size of 
its error bars. For the purpose of estimating the size of 
the background to be subtracted, the level of agreement 
is more than adequate. The earlier discussion for the 
12 C R p - wave distribution applies here as well; the data 
are consistent with a negligible background contribution. 
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C. Self-consistency checks 

It is of interest to investigate the validity and self- 
consistency of the analysis followed above. Our method 
is to check whether the fit coefficient values, used to de- 
termine the background contribution to be subtracted, 
are able to reproduce the cos^* + data distributions for a 
number of different bins. For example, if R p - W ave 
was not a sufficiently accurate measure of the shape of the 
data's cos6** distributions, the A, B, C coefficient values 
determined from fits to them would not provide good 
descriptions of the data, and the validity of the back- 
ground subtraction just performed would be in doubt. 
The 800-960 MeV photon energy region was chosen for 
this check, as it has the best statistical precision, and 
so any systematic deviations from the data may be most 
easily discerned. 

Fig. 1191 shows the distributions for four m ffT bins cen- 
tered at approximately 552, 602, 646 and 695 MeV/c 2 . 
As discussed earlier, each invariant mass bin is 90 
MeV/c 2 wide, and so there is partial overlap between 
adjoining bins. The long-dashed line in each panel is the 
best fit to the data with the A + Bcos9* + + Ccos 2 0* + 
function. The short-dashed line is a result of the sum 
of the p° L and background contributions, as discussed in 
the previous section. The agreement of the two curves is 
a direct measure of the reliability of the identification of 
the foreground and background components in the data. 
For all three nuclei, the curves agree quite well through- 
out the m TT range. While it is to be expected that the 



FIG. 19: Helicity angle distributions, with opening angle 
and missing mass cuts appropriate for each nucleus applied, 
for four invariant mass bins and 800 < E~j < 960 MeV. The 
legend at right indicates the mean invariant mass value of the 
data in each bin, averaged over the three target nuclei. The 
long-dashed curves are the fits to the data used to extract 
the R p -wave distributions. The short-dashed curves use the 
coefficient values determined from the fits to the Rp-wave 
values presented in Figs. 1151 1161 and 1171 The consistency of 
the two curves, within data error bars, indicates the validity 
of the p-wave shape ratio analysis. 



long-dashed curves will be better descriptions of the data, 
as they have been directly fit to them, the deviations 
of the short-dashed curves from the data are compara- 
ble to the size of the data error bars. The 602 and 646 
MeV/c 2 bin data have the smallest statistical error, and 
the two curves are in the best agreement. The 552 and 
695 MeV/c 2 bins have significantly fewer events, and the 
discrepancy between the curves is larger, as expected. 
The overall conclusion is that there is no systematic de- 
viation of the short-dashed curves from the data which is 
outside the statistical error of the data and analysis, and 
so the reliability of the background subtraction method 
is validated. The overall self-consistency achieved for the 
12 C data is quite impressive. This is a direct result of the 
very effective suppression of the background processes by 
the two cuts. 
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VI. COMPARISON OF THE EXTRACTED 
IN-MEDIUM pi INVARIANT MASS 
DISTRIBUTIONS WITH MODEL 
CALCULATIONS 

Before our extracted invariant mass distributions can 
be compared to any model, a number of effects must be 
taken into account in its prediction. Because of the sub- 
threshold nature of the experiment, the role of the lim- 
ited kinematic phase-space is significant, and must be 
accounted for. In addition, although every effort has 
been made to quantify and remove non-p^ background 
contamination from the data, the effect of the TAGX 
spectrometer acceptance and detection thresholds still 
remain. Similarly, the model prediction should be sub- 
jected to the same O^n and m m i SS cuts as the experimen- 
tal data. The best way to take these effects into account 
is to embed the model calculation within a MC simula- 
tion of the experiment, and analyze the simulated events 
in the same manner as the data. Only by following a pro- 
cedure of this manner can reliable comparisons be made 
to the data. 

In this section, we will first present a comparison of our 
data to two kinematic models utilizing the free p° line- 
shape obtained from PWA analysis of the e + e _ — * tt + it~ 
reaction |23j . This free line-shape provides a superior 
description of the low mass p° tail, significantly better 
than the standard parameterization of Ref. |2(|. This 
will investigate the roles of kinematic phase-space and 
Fermi momentum and will provide the first indication 
of the reaction mechanism responsible for the significant 
longitudinal p° polarization observed here. Then, we will 
compare our data to three phenomenological models of 
the in-medium p® L . We believe this will provide the most 
rigorous comparison to date between model and obser- 
vation of the in-medium p° characteristics. In all cases, 
the simulated events were tracked through the simulated 
TAGX spectrometer and have all acceptance and thresh- 
old effects applied. The simulated events were then an- 
alyzed in exactly the same manner as the data, and the 
quadratic fit coefficient to the cos 2 6>* + distribution ex- 
tracted, just as for the data. 



A. Comparison with the jNf — > pi N kinematic 
model 

In Fig. 1201 the extracted invariant mass distributions 
for all three nuclei are shown together with the quasi- 
free pi MC simulation discussed earlier [long-dash lines] . 
Overall, we see that the simulation provides a good de- 
scription of the data for all three nuclei and all tagged 
photon energy bins. This is a convincing statement that 
indeed the p-wave signatures isolated are those of sub- 
threshold p L production. The probability that uncor- 
rected 7T + 7r _ pairs from unrelated production processes 
could result in an invariant mass distribution similar to 
the partial-wave analysis p° and respond in the same 



manner to the applied kinematic cuts, is simply too re- 
mote to be credible. 

Upon closer inspection, we observe deviations between 
the data and quasi- free p° model distributions. The 
agreement between data and simulation is best for 12 C. 
For the 600-800 MeV E 7 bin, the agreement between 
the 12 C MC simulations and the data is excellent for 
the whole range of the invariant mass values. For the 
800-960 MeV bin, the width of the simulated distribu- 
tion is the same as that of the data, but it is peaked at 
slightly higher mass. In the case of the 960-1120 MeV bin, 
the quasi-free p° L distribution is slightly broader than the 
data and peaked at a higher invariant mass. It appears 
that the trend barely discernible in the 800-960 MeV bin, 
has become more pronounced in this energy region. 

The same broad conclusions are also true for the com- 
parison between MC and data for 3 He. Here, the widths 
of the data and MC distributions are nearly the same, 
but the data are consistently peaked at about 20 MeV/c 2 
lower mass. Surprisingly, the deviations between data 
and simulation are greatest for 2 H. In both the 600-800 
and 800-960 MeV bins, the data distributions are nar- 
rower than those of the simulations, and there is a pro- 
gressive shift in the centroids of the data distributions 
away from the simulations as the photon energy is in- 
creased. 

In conclusion, while there is qualitative agreement be- 
tween the data and the calculated invariant mass distri- 
butions, the deviations merit further study. The higher 
centroids and broader widths of the quasi-free MC distri- 
butions may indicate more phase-space available to the 
assumed production and decay reaction than the data 
justify. 



B. Comparison with the 7iV F -> iV*(1520) -> p%N 
kinematic model 

Fig. |20| also compares the data distributions to the 
quasi-free TV* (1520) — > p° L N mechanism MC presented 
earlier [short-dash lines] . The assumed free p° line-shape 
is the same, but the effects of kinematic phase-space and 
Fermi momentum arc different, due to the more con- 
strained reaction mechanism. For the 2 H 600-800 MeV 
energy bin, the agreement between data and the simula- 
tion is much improved from the quasi-free p a L MC and it is 
now quite impressive. The different phase-space available 
to the resulting p° L has narrowed the calculated distribu- 
tion and the overall agreement extends throughout the 
TO W7r range. It is the strongest evidence yet that TV* (1520) 
decay is the dominant p° L production mechanism in this 
energy range for this nucleus. 

The improved agreement between the iV*(1520)-based 
production mechanism and the 2 H data carries over to 
the 800-960 MeV energy range. In comparison to the p° L 
MC, the simulated distribution has a more similar width 
to the data, although the centroids disagree. For the 960- 
1120 MeV energy bin, the data and TV* (1520) MC distri- 
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butions appear shifted by approximately 40 MeV/c 2 from 
each other, with the data being lower than the model dis- 
tribution. The precise reason for this is not clear but it 
may be related to the quasi-free TV* (1520) assumption, 
admixture of coherently produced p° and/or admixture 
of quasi-free p\ production. We conjecture that this ad- 
mixture may also be responsible for the small bump near 
m W7r = 800 MeV/c 2 in this energy bin. 

Now moving on to the 3 He distributions, we see that 
the iV*(1520) production mechanism gives a much better 
description of the 800-960 MeV data than the quasi-free 
p\ simulation. As expected from the experience with 
2 H, the iV*(1520)-based distributions are both narrower 
and shifted to lower values compared to the p\ quasi- 
free simulation. For the 960-1120 MeV energy region, the 
3 He data are also in better agreement with the A* (1520) 
process, but the level of agreement is less impressive than 
for the lower energy bin. 



Finally, for 12 C the level of agreement with the 
iV*(1520)-based production mechanism is worse than 
with the quasi-free p\ mechanism for the 600-800 and 
800-960 MeV bins. For the 960-1120 MeV bin, the cen- 
troid of the distribution is given better by the TV* (1520) 
simulation, but the width is given better by the quasi-free 
p° L simulation. While the two lighter nuclei are broadly 
consistent with dominance of the TV* (1520) channel, for 
12 C, only in the 960-1120 MeV energy regime can some 
evidence for the iV*(1520) channel be found. 



C. Conclusions from comparison to kinematic 
models 

The general agreement of the experimental p Q L invari- 
ant mass distributions with the expectations from the 
PWA-based parameterization verify that the analysis has 



21 



indeed identified pions originating from p° decay. The 
small Fermi momentum distribution of a single nucleon 
in 2 H allows sub-threshold production of the p°, while 
the small nuclear density is expected to result in min- 
imal, if any, medium modifications. As such, 2 H is an 
ideal nucleus to use as proof of principle of the analysis 
that indeed the p° has been identified and been isolated 
in the event sample. This analysis, then, should lay to 
rest any question whether the p° has been indeed identi- 
fied at sub-threshold energies. 

The 2 H and 3 He data invariant mass distributions favor 
p\ production via the quasi-free TV* (1520) — > p\ TV mech- 
anism. The two quasi-free mechanisms investigated here 
have invariant mass distributions that are significantly 
different from each other, allowing these comparisons to 
be meaningful. The comparison with the kinematic mod- 
els has thus given a strong experimental indication of the 
likely p\ production mechanism at these photon energies. 
Regarding the lack of evidence in the 12 C data to sup- 
port even a modest content of the latter mechanism, we 
note that the role the TV* (1520) plays in nuclei, and its 
in-medium mass and width, is an ongoing debate. Our 
simulation made use of free TV*(1520) parameters from 
Ref. |2(| . While there is good evidence to suggest that 
the TV* (1520) is significantly broadened in complex nu- 
clei, more recent results from photo-absorption studies on 
carbon and more massive nuclei shed new light on this 
subject [35| . Even though the photon energies did not 
exceed 800 MeV in that work, its primary subject was 
the TV* (1520) excitation and its possible medium mod- 
ifications, including broadening, in nuclei. The lack of 
experimental signatures of the resonance in total photo- 
absorption cross sections in nuclei like carbon, thus re- 
mains a not-understood phenomenon. 

Finally, a detailed comparison of the data with the sim- 
ulation results indicate deviations from the simple kine- 
matic models considered here, which are likely due to 
higher TV* and A resonances and/or reaction mechanism 
admixtures. The comparison of the data with the phe- 
nomenological models to follow may shed further light on 
these discrepancies. Of the three models presented, two 
are based on established hadronic interactions that do 
not involve quark or QCD elements. While both models 
share a similar philosophy and share some common el- 
ements, there are substantial differences between them. 
The third model is also phenomenological in nature, but 
it is based on quark degrees of freedom and its founda- 
tions are thus completely different from the other two. 



D. Comparison with the Rapp - Chanfray - 
Wambach (RCW) Model 

This model p| has two main contributions to the in- 
medium spectral shape of the p. First, the interactions 
of the pions with the surrounding nucleons and A's accu- 
mulate substantial strength in the p spectral function for 
the lower regions of the invariant mass. This renormal- 



ization of the pion propagation in a 7rTVA gas was first 
proposed in j^o, Ej} ■ Second, the contributions to the p 
spectral shape from in-medium p-baryon scattering are 
evaluated. This model reproduces the experimental data 
on p-wave tttt scattering in free space as well as the pion 
electromagnetic form factor in the time-like region. 

The first element of the model accounts for medium 
modifications of the p by considering the interaction of 
the intermediate two pion states, in other words, the p is 
"dressed" by the two pion intermediate state. The single 
pion self energy in these intermediate states is evaluated 
within the particle-hole excitations at finite temperature. 
Specifically, the pions interact with the surrounding nu- 
cleons and with the thermally excited A states through 
excitations of the type TVTV -1 , A TV -1 , AA" 1 and TVA -1 . 
Evaluation of these type of interactions lead to substan- 
tial broadening of the p spectral function. The peak is 
shifted to slightly higher values for the invariant mass, 
however, significant strength is added to the low mass 
region below 600 MeV/c 2 . 

The second element of the model treats the scattering 
of the p in-medium. The sizable strength of the pNN and 
pNA coupling constants, and the large branching ratios 
of TV(1720) and A(1905) decays to pN final states, led 
the authors to consider pTV(1720)TV _1 and pA(1905)TV" 1 
particle-hole-like states called "rhosobars". Such states 
have been extended to include the A(1232) resonance. 
The calculated contributions of the overlapping states of 
such "rhosobars" result in an appreciable enhancement of 
the invariant mass distribution below 600 MeV/c 2 while, 
at the same time, leading to a depletion of the p invariant 
mass peak. 

These two main sources of p medium modifications lead 
to an enhancement of the two-pion invariant mass spec- 
trum below the region accounted for by free p and u) 
mesons and they appear to account very well for the ex- 
perimental data on di-lepton production in high energy 
heavy ion collisions @ and [2a] • As such, the model pro- 
vides an alternate explanation of the di-lepton spectra 
to that of the chiral phase transition [39l Eoj . Although 
the model is particularly suited to conditions of nuclear 
matter at high densities and temperatures (heavy ion col- 
lisions at high energies), it can also find application to 
nuclear matter at lower densities such as those found in 
the nuclear core. It is a non-relativistic model and this 
is the main difference from the model of Post - Leupold 
- Mosel P. 

Details on how the RCW-model was implemented in 
the TAGX simulation are discussed in App. IA II The 
RCW-model simulation is compared to the data in Fig. 
1211 In comparison to the two kinematic models in Fig. 
1201 the RCW model provides a good description of the 
data over the full energy range and for all three nuclei. 
Its description of the 12 C data is similar to the quasi- free 
p\ kinematic model while at the same time its descrip- 
tion of the 3 He data is similar to the TV*(1520) kinematic 
model. For the 2 H data, it does better than the p° L quasi- 
free model, but still not as good as the TV* (1520) model. 
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This could be due to the limited number of spectral func- 
tions provided (see App. IA II) . which are insufficient to 
describe the high density regions of the deuteron. As 
with the kinematic models, the agreement with the data 
is worst for the 960-1120 MeV bin. The model curves 
give broader distributions than the data justify and this 
is perhaps indicative of the lack of the A* (1520) in the 
model. Although A's and A*'s are included, the latter 
are of higher mass than the iV*(1520) and play little, if 
any, role at our relatively low photon energies. 

E. Comparison with the Post - Leupold - Mosel 
(PLM) Model 

This model is made in the low-density approximation, 
where the in- medium self energy of the p is completely de- 
termined by the pN forward scattering amplitude. This 
is a critical element of the model and the pN forward 
scattering amplitudes are extracted from irN — ► irirN 
processes ^^J. The model is relativistic and it thus 
avoids the problems of various components of the the- 
ory being evaluated in different frames of reference, as in 
other phenomenological models. In addition, the prob- 
lem of the non-coupling of longitudinal p states to p-wave 
resonances in non-relativistic models is resolved in this 
model. 

An essential ingredient of this model is the connec- 
tion via unitarity between the matrix element for the 
decay of a resonance R into a nucleon and a p of a given 
polarization \M^^ p \ 2 to the imaginary part of the for- 
ward scattering amplitude T T / L : Im T T / L ~ \M^\. 
All the baryon resonances with sizable couplings to the 
pN channel have been included with relevant parame- 
ters from |45|. For transversely polarized p mesons, the 
A''* (1720) resonance plays a dominant role. In the lon- 
gitudinal channel, the A* (1520) is the main contributor. 
For some of the resonances, there is considerable vari- 
ation among different parameterizations of the coupling 
strength to the pN decay channel. 

This model concludes that the A* (1520) plays a large 
role on the propagation of p mesons in nuclei. While 
in the longitudinal channel the differences between rela- 
tivistic and non-relativistic calculations are small, in the 
transverse channel and at large momenta there are sig- 
nificant differences and this model leads to pronounced 
broadening of the p invariant mass peak. 

Details on how the PLM-model was implemented in 
the TAGX simulation are discussed in App. IA 21 The 
PLM model authors felt that the model could only be 
applied reliably to 12 C, and so a comparison will only 
be made to the data for that nucleus, in Fig. For 
the 600-800 MeV tagged photon energy bin, the PLM 
model provides a good description of the width of the 
data distribution, but the centroid is shifted to slightly 
lower values. No other model has under-estimated 
the centroid of this data distribution. The agreement be- 
tween the PLM model and the 800-960 MeV data is im- 



proved, but is not nearly as good as that provided by the 
RCW model. The agreement with the E 1 = 960 - 1120 
MeV data is not as good as for the mid-energy region, 
with the model distribution being broader and shifted to 
higher invariant mass values. 

Nuclear absorption of one or both pions may distort 
the spectral shape due to different absorption rates in 
different density regions of the nucleus, while medium 
modifications also depend on nuclear density. The effect 
of pion absorption in the medium was considered via the 
method described in section ITlI A 41 and is shown via the 
dotted curve in the figure. The simulation indicates that 
after the m m i SS and 9^ cuts are applied, the effect of 
pion absorption in 12 C is to remove events in an equal 
proportion from all parts of the m^n distribution, re- 
sulting in an invariant mass distribution which is almost 
indistinguishable from the simulation in which pion ab- 
sorption is not considered. 



F. Comparison with the Saito - Tsushima - 
Thomas (STT) Model 

The foundations of this model lie in the quark-meson 
coupling (QMC-II) version of the original QMC model 
of Ref. O E3 utilizing the MIT quark-bag model. 
Here, quarks in non-overlapping nucleon bags interact 
self-consistently with both scalar and vector mesons in 
the mean-field approximation. The vector mesons are 
themselves described by meson (quark) bags, and the 
bag parameters are fixed to reproduce the free nucleon 
mass and radius. Since the vector mesons are described 
by quark bags as well, additional parameters have been 
introduced to fit their free (vacuum) masses. The in- 
medium masses are then given in terms of the mean-field 
value of the a meson at that density, which is given by a 
parameterized (phenomenological) function. 

It is worth noting that no pion or p° interactions with 
the nucleons in the nuclear medium, A and N* excita- 
tions included, are incorporated in this model. The p° 
interacts with the scalar field of the nucleus and the two 
pions emerge with their original four momenta. Thus, 
this model provides a description of the modification of 
the p° mass due to nuclear binding effects only. 

Details on how the STT model was implemented in 
the TAGX simulation are discussed in App. IA 31 The 
simulation is compared to the data in Fig. [231 There is 
a slight shift between the 2 H and 3 He data and model 
distributions, but overall the agreement is impressive, 
especially given the model's lack of any in-medium p° 
width or shape modification. Ref. |43j was the only 
work to provide specific predictions for our experiment; 
their predicted result for 3 He was a 40 MeV reduction in 
the mass of the p°. The mean of the 800-960 MeV data 
distribution is lower than the predicted m* value by ap- 
proximately 20 MeV. The STT model is the only one 
to provide a good description of the 3 He data distribu- 
tion at 960-1120 MeV. As the model has no sensitivity to 
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FIG. 21: Invariant mass distributions of the p° L strength extracted from the data [symbols] compared to quasi-free simulations 
[short dash-long dash lines] incorporating the Rapp-Chanfray-Wambach model Q in-medium p° L spectral function. In the 
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both the statistical precision of the data, as well as the uncertainty in the background subtraction. 



o 
. f— i 
■+-> 

• r— I 
■+-> 

w 
Q 

H 



600-800 



1.2 
1.0 
0.8 H 
0.6 
N 0.4 
0.2- 



0.0 



■M 

I 1 



i 1 



400 500 600 700 



800-960 



1.2 
1.0 
0.8 
0.6 
0.4 
- 0.2- 



0.0 



AT— ! 



jj 



X I 



500 



600 



700 



1.2 
1.0 

- 0.8- 
0.6 
0.4 

- 0.2- 
0.0 



960-1120 



800 



600 



-at. 



in 



700 



800 



900 



m Invariant Mass (MeV) 



FIG. 22: Invariant mass distributions of the p° L strength extracted from the data [symbols] compared to quasi-free simulations 
[dash-dot lines] incorporating the Post-Leupold-Mosel model |f| in-medium p° L spectral function. The dotted lines indicate the 
result when 7r absorption is added to the simulation. The data error bars reflect both the statistical precision of the data, as 
well as the uncertainty in the background subtraction. 



24 



spectral shape modifications, this agreement would imply 
that the p° L line shape is largely unaffected by the nuclear 
medium, a fact consistent with the good agreement with 
our PWA-based kinematic models. 

For 12 C, the agreement between model and data is also 
quite good. At 800-960 MeV, the theoretical curve has 
shifted to slightly lower values compared to the data, and 
has a narrower distribution. This is the only case where 
a model distribution is sharper than the data indicate. 
This trend becomes more pronounced in the 960-1120 
MeV range and, at the same time, the distribution has 
become significantly broader than the data. This behav- 
ior is unlike all other theoretical or calculated responses 
and it is unique in its energy dependence. It is inter- 
esting to note that, in comparison to the 960-1120 MeV 
data distribution, the STT and RCW-PLM models ex- 
hibit shifts in opposite directions. A model which incor- 
porates both the nuclear binding effect of STT as well 
as the spectral function modifications of RCW-PLM has 
the potential to provide the best description of the 12 C 
data. 



G. Conclusions from comparison to 
phenomenological models 

The first model compared to the data, RCW, provides 
a reasonable description of the data from all three nuclei 
over all energy bins. Despite the fact that the V*(1520) 
is not explicitly included in the RCW model, the invari- 
ant mass distributions from this simulation exhibit some 
of the characteristics of the quasi-free p° and quasi-free 
V*(1520) kinematic simulations, presumably a reflection 
of the "rhosobar" nature of the model. 

The second model, PLM, is more sophisticated, as it is 
fully relativistic. Where they can be directly compared, 
the differences between the RCW and PLM model simu- 
lations are small. While the RCW-simulated distribution 
for 12 C tracks the changing shape of the data distribution 
with energy, the PLM model predicts a larger energy de- 
pendence than the data portray. The better agreement 
between the PLM model and the data at 800 < E 1 < 960 
MeV is due, in part, to the fact that it is bracketed by 
two energy ranges where the simulated distributions are 
respectively at too low and too high mass compared to 
the data. Given the inputs and complexity of both mod- 
els, the identification of possible factors that allow one 
model to give a better description than the other is not an 
easy task. However, the clear superiority of the quasi-free 
p\ kinematic simulation over the V*(1520) for 12 C may 
provide a helpful indication of why the RCW model does 
better. The RCW model does not include the JV*(1520), 
while in the PLM longitudinal channel model it has sig- 
nificant importance. Since the higher mass resonances in 
both models are largely suppressed due to phase-space in 
this experiment, we conjecture that the PLM model may 
over-estimate the role the iV*(1520) plays in 12 C. 

Given its simplicity, the STT model does a surpris- 



ingly good job of describing the data. Because the free p° 
width is assumed in our simulation, the good agreement 
of the STT-simulated distribution with the 3 He data, for 
example, raises the possibility that what is observed for 
this nucleus is the combination of a free p° spectral func- 
tion and a central mass shift. For 12 C, the STT model 
under-predicts the p° mass, and it might seem reasonable 
that the combination of the STT and RCW-PLM mod- 
els will together provide a good description of the data, 
i.e. a central mass shift in addition to a spectral function 
modification. 



VII. COMPARISON WITH OUR PREVIOUS 
RESULTS 

This work reports new results on 2 H and 12 C and a re- 
analysis of our 3 He data in light of the helicity anal ysis 
reported in Ref. 01- O ur earlier work on 3 He [ToL IllL 
Il2ul.3| reported significant mass modifications and, in the 
process, it proved to be controversial. It is of interest then 
to compare the older results with the analysis reported 
here and investigate their differences. 

The 3 He analysis in Refs. [H IH El 

was based on 

the simultaneous fitting of a number of quasi-free sim- 
ulations resulting in n + Tr~ production to data for five 
different kinematic observables. Ref. ^3 was f° r the 
tagged photon energy range of 380-700 MeV, while Refs. 
[TfllH reported results for 800-1120 MeV. In compar- 
ison to this work, these analyses used relatively loose 
cuts to enhance the p° contribution relative to the back- 
ground channels. For example, the analysis of Ref. ^3 
was based on a simultaneous fit to the observables from 
the entire data sample (no physics cuts) as well as to a 
set of observables to which a set of p° enhancement cuts 
(70° < 6U < 180° and 2700 < m misa < 3050 MeV/c 2 ) 
were applied. Because nearly the entire data sample was 
used in the analysis, a mass modification result could 
even be extracted for the E 1 < 700 MeV range, while this 
was not possible for the more restrictive helicity analysis 
reported here. On the other hand, because the conclu- 
sions were explicitly based on the fit of a series of MC 
simulations to the data, the results of Ref. 0, 0, 
are more model-dependent than the helicity analysis re- 
ported here. So, the two independent analyses have their 
respective strengths and weaknesses, and taken together 
they can provide a more complete view of the properties 
of the in-medium p° meson. 

Refs. mmm 

restricted their attention to the 'best 
fit' p° effective mass value, after accounting for back- 
ground process contributions, free-mass p° production, 
and kinematic phase-space effects. Since the emphasis 
was only on the extraction of the most-probable reduced- 
mass value, a Gaussian line shape was taken for the p , 
using the free width. This choice of line shape has been 
criticized, but the use of a Lorentzian line shape was 
investigated in Ref. ^2 an d yielded almost identical re- 
sults. Even though these two line shapes have quite dif- 
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FIG. 23: Invariant mass distributions of the p° L strength extracted from the data [symbols] compared to quasi-free simulations 
[dashed lines] incorporating the Saito-Tsushima-Thomas model |4*ll in-medium p° L mass. In the panels for 12 C, the dotted lines 
indicate the result when 7r absorption is added to the simulation. The data error bars reflect both the statistical precision of 
the data, as well as the uncertainty in the background subtraction. 



ferent low and high m TI tails in unresctricted (free) phase 
space, when the TAGX thresholds and E 1 imposed phase 
space are taken into consideration, the resulting distribu- 
tions are nearly identical. The conclusion of these studies 
was that the "effective" p° mass in 3 He was reduced sub- 
stantially from its free value. 

Ref. [3 found an effective mass of 490 ± 40 MeV/c 2 
for 460 < Ery < 700 MeV. However, no further compar- 
isons to the present analysis can be pursued. For the 
higher energy 3 He data, Ref. ^(| found the largest im- 
provement in the relative x 2 fit to the kinematic distri- 
butions at 800-880 MeV, with the next largest improve- 
ment for the 880-960 MeV bin. Ref. [nj, which in- 
cluded several improvements, concluded that these two 
energy regimes are the only ones where mass modification 
could be reliably extracted, with effective mass values of 
m* = 642 ± 40 MeV/c 2 for 800 < E 1 < 880 MeV and 
77i* = 669 ± 32 MeV/c 2 for 880 < E 1 < 960 MeV. The 



960-1040 MeV data were marginally supportive of an ef- 
fective mass value of m* = 682 ± 56 MeV/c 2 , whereas 
the data in the 1040-1120 MeV region were consistent 
with no p° mass modification. A change in the width of 
the p° was briefly investigated in Ref. and it was 
concluded that "these exploratory fits verify the prefer- 
ence for a reduced p° mass, but are inconclusive, within 
the sensitivity of the data, as to whether a width modifi- 
cation is supported in addition. " Finally, in Ref. |48[ we 
suggested that the observed effective p° mass could either 
be due to the proximity of the p° to the struck nucleon 
and the shorter p° mean decay length at lower energy, or 
due to the collective excitation of a iV*(1520)iV~ 1 state. 

An independent analysis of the 3 He data was published 
m Ref. 12]. Cuts were placed to identify the exclusive 
p° production (m m i ss ~ rn3 He ) and the quasi- free p° 
production regions (m m i SS < (m3 He + m,) for further 
analysis. The p° line-shape was taken as a Breit-Wigner 
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form and non-p° background was modeled via a phase- 
space mechanism. The interferences between free-mass 
p°, reduced-mass p° and non-p" contributions were taken 
into account via the Soding 49] model. For the exclusive 
production region, the extracted invariant mass distri- 
bution was consistent with free p° parameters, presum- 
ably these events originated near the periphery of the nu- 
cleus. The analysis of the quasi-free production region, 
expected to correspond to production deeper within the 
nucleus, yielded a line-shape significantly at odds with 
the free p° parameters, implying a reduced-mass value of 
m* = 655 ±8 MeV/c 2 . 

For comparison with the present work, we restrict our 
attention to 3 He. The good agreement of the free PWA 
line-shape-based simulations with the data (e.g. Fig. I20|) 
support the conclusion that the p° width has not been 
substantially changed from its free value. This is consis- 
tent with the theoretical expectations of the RCW and 
PLM models, where the dominant width modification is 
restricted to transversely polarized p 1 ^- Comparison of 
the PWA line-shape-based simulations with the 960-1120 
MeV data support a mass shift. Similarly, comparison of 
the RCW model with the data (Fig. EH), indicate 
that the data support a slightly lower mass value than 
the model predicts. This conclusion is also reinforced by 
the comparison of the STT model 43] with the data (Fig. 
12 3|) , which support a reduction of the p° mass while using 
the free width. The STT model predicts [3 a 40 MeV re- 
duction in the p° mass in 3 He, and essentially agrees with 
the 960-1120 MeV data in Fig. EH confirming to* = 730 
MeV/c 2 . For the 800-960 MeV bin, the STT model dis- 
tribution is ~ 20 MeV high compared to the data, giving 
a reduced mass value no higher than m p o ~ 710 MeV/c 2 . 
Although these values are obtained with reference to the 
STT model, they are consistent with our expectations 
from the kinematic and RCW models, and so are listed 
in Table IIIII as representative 'effective masses' for this 
work. It is difficult to place an error bar on these masses 
within the context of the STT model, and so we refrain 
from doing so. 

TABLE III: m* results for 3 He from this work and from Ref. 
[ll|. For this work, the masses shown are determined via 
comparison to the STT model |43| simulation. Ref. [Tl| pre- 
sented the data contributions from 'free mass' and 'reduced 
mass' p° components separately, and so they must be com- 
bined in a yield-weighted manner before comparison to the 
results of this work. 





STT model m* 


Yield-weighted m* 




This Work 


Ref. [11] 


800-960 MeV 


700-710 MeV 


672 ± 31 MeV 


960-1120 MeV 


730 MeV 


743 ± 17 MeV 



To compare the old and present analyses, we must not 
only use the same tagged photon energy bin limits, but 
we must also keep in mind that the old analysis separated 
the free-mass and reduced-mass p° components, while 



the present results are for the total observed p\ yield. 
At 800-960 MeV, free-mass p° production is very small, 
consistent with the 971 MeV threshold for p° produc- 
tion one a below nominal mass on the proton, and then 
its contribution grows rapidly thereafter. Yield-weighted 
averages of the Ref. [ll| results, including the observed 
free-mass p° component, are listed in Table ITTT1 

Given their substantially different analysis techniques, 
the two sets of results in Table lTTll are in broad agreement. 
The to* values from the present work are slightly higher 
than our earlier results from Ref. 11], but agree within 1 
a. While the results of the present work are intrinsically 
less model-dependent, they rely on tight cuts upon the 
data, and so we were unable to extract any result for the 
600-800 MeV energy bin, unlike our earlier works. This 
work verifies the conclusions from Refs. |ToL fTTL f]~3 | that 
there is significant medium modification observed on the 
3 He target as a result of the sub-threshold production 
technique. A second conclusion regards the issue of p° 
width modification. While our old analysis was unable 
to make any statement on this issue, the agreement of 
the PWA-based simulations with the data distributions 
appear to rule out any substantial p° L width modification, 
in accordance with model expectations. 



VIII. DISCUSSION AND CONCLUSIONS 

Employing the unique signature of longitudinally po- 
larized p°r in the cos6** + distributions established in Ref. 
|14| , this work has identified and isolated p L events from 
tagged photo-production on 2 H, 3 He and 12 C in the 
E 1 = 600 — 1120 MeV energy region. The analysis led to 
the extraction of the invariant mass distributions of the 
p\ for all three nuclei over three bins of tagged photon 
energy. For 2 H and 12 C, this is the first time such re- 
sults are reported. The 3 He data were re-analyzed in the 
same manner and the results compared with our previ- 
ously published results in Refs. [lOt fill LL2| . 

The effects of the two cuts used in the helicity analy- 
sis have been studied by comparison to a series of MC 
simulations with the same cuts applied. It has been 
shown that these cuts cannot create artificial longitu- 
dinal helicity signatures mimicking a p Q L contribution to 
the data. A p-wave-like component from non-p^ resid- 
ual background was identified with the assistance of a 
p-wave shape ratio and then subtracted, yielding the ex- 
perimental in-medium p° L [m^ , E 7 ) distributions. This 
background component was found in most cases to be 
quite small, and its subtraction yielded distributions 
nearly indistinguishable from the unsubtracted distribu- 
tions. Thus, the systematic uncertainty in the extracted 
p\ distributions due to the background subtraction is 
small. 

The 600-960 MeV 2 H target data favor production 
via the -fN F -> V*(1520) -> p° L N' production mecha- 
nism. For the 960 < E 1 < 1120 MeV bin, neither of the 
kinematic production mechanism models account for the 
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observed distribution. The agreement of the data with 
models based on the PWA-line-shape for the free p° is a 
conclusive statement of the identification of the p\ events 
in the data sample, and of the reliability of the analysis. 

The 3 He distributions are also consistent with produc- 
tion via the A* (1520) mechanism, over the full observed 
photon energy range of 800-1120 MeV. Two phenomeno- 
logical models 0, are quite successful in describing 
the data. Comparison with the STT model simulation, 
using the free p\ width and line-shape indicates that the 
data provide no evidence of an in-medium p° L distribu- 
tion that is broader than the free distribution. This is 
consistent with model expectations, which predict that 
the in-medium p\ distribution will be significantly wider 
than the p\ distribution. The STT model 0] predicts a 
modified mass of approximately 730 MeV/c 2 in 3 He, and 
it is in excellent agreement with the 960-1120 MeV data. 
For the 800-960 MeV bin, the data support a somewhat 
lower mass, to*„ = 700 — 710 MeV/c 2 . These results are 

slightly higher than the older analyses of the same data 
in Ref. |ToL ITTl Il2| , but they confirm the observation of 
medium mass modifications in the 3 He case. 

The 12 C distributions are consistent with quasi- free p\ 
(PWA) production over all three tagged photon energy 
bins. The data are in excellent agreement with the RCW 
model 3 for the two lower energies, but RCW predicts a 
wider distribution in the 960-1120 region than the data 
support. The PLM model [f| does well overall, but re- 
sults in a larger variation in the distribution with 
energy than the data show. The STT model [43 pre- 
dicts distributions narrower than the 12 C data for the two 
lower energy bins and it is the only one of the three phe- 
nomenological models that predicts a lower mass value 
than the data exhibit. This indicates that a model which 
incorporates both the nuclear binding effect of STT as 
well as the spectral function modifications of RCW-PLM 
has the potential to provide the best description of the 
12 C data. 

This very line of argument has already been raised 
by Brown and Rho |50j in the context of recent re- 
sults from RHIC. The Brown-Rho scaling of medium- 
dependent masses m*/m p ~ f^/fn is combined with the 
spectral function modification of RCW in a unified pic- 

[55] The data range of |t| values given in Ref. Il4l is incorrect, the 
actual range displayed by the data is 0.2 to 1.1 (GeV/c) 2 . 



ture. In this case, the p° mass in the phenomenological 
Lagrangian used by RCW is replaced with what is es- 
sentially the Brown-Rho scaling mass m*. For the tem- 
peratures and densities appropriate to this work, the two 
mechanisms push the p°-mass in the same direction, and 
so result in a lower p° mass than either model alone pre- 
dict, nearly doubling the density effect [5l| . 

The physical mechanism for the dramatic helicity flip 
remains an open issue. These data do not support either 
transversely polarized pj, or unpolarized p° production. 
Either would have manifested as a cr°-like background 
associated with some strength at the central region of 
the cos0* + distribution. In Ref. 14], one explanation of- 
fered for the observed longitudinally polarized p° L was the 
large \t\ in that work associated with p° production [55] 
Rapidly increasing helicity-flip amplitudes had been ob- 
served in 0] for values of \t\ larger than 0.4 GeV/c. An 
alternate explanation is that the sub-threshold reaction 
mechanism may be responsible for the observed helicity 
flip. 

Finally, the production mechanism seems to depend 
on the nucleus in question. Even though A* (1520) exci- 
tation in 2 H is both broader and quenched with respect 
to the free proton, it constitutes a well-identified con- 
tribution to the second resonance group in pion photo- 
production studies [35) 1 (and references therein). For 
the Ej energies of this work, one expects significant 
A''* (1520) excitation. The strong decay channel into p° 
|26| makes this a favorable production mechanism of the 
low invariant mass part of the p° . The data are consistent 
with this expectation. 
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APPENDIX A: THE IMPLEMENTATION OF 
THE PHENOMENOLOGICAL MODELS IN THE 
TAGX SIMULATION 

1. The Rapp-Chanfray- Wambach model 

In-medium p° spectral functions were obtained 
for cold nuclear matter at three matter densities, p' B = 
Pb/Puuc = 0.5, 0.7 and 1.0, where p nuc indicates the 
standard nuclear density of 0.155 fm -3 . Longitudinal, 
transverse and spin-averaged functions were provided, 
but as the data strongly support the production of p® L , 
only the longitudinal spectral function was used here. 
The in-medium p° propagator for each density is given 
in matrix form as functions of invariant mass and three 
momentum q, where q is the momentum of the p° with 
respect to the nuclear medium, 

Im = Im Dp(p' B ,m^,q). 

For each simulated event, the tagged photon energy, p° L 
mass m n7r , and reaction location within the nucleus were 
randomly chosen and an event generated via the quasi- 
free mechanism jNp —> p® L N , where Np is the partic- 
ipating struck proton with initial Fermi momentum pp 
and the remainder of the nucleus is a spectator. It should 
not be necessary to consider any other kinematic chan- 
nel, such as N* — > p°N, as these are already taken into 
account via the in-medium spectral function. To do oth- 
erwise would lead to "double counting" of resonance ef- 
fects. The Fermi momentum distributions of [28l.l29l.l3fH 
were used. 

The nuclear density at the production location was cal- 
culated in the following manner: 

12 C: The electric charge distribution data of Ref. |42] | 
was parameterized as 

p{r) = a (l + a 1 r 2 )e a2r2 p n uc 

where ao = 1.2550, a\ = 0.53765 fm -2 , = 
—0.39486 fm -2 , and r is in fm. Note that the stan- 
dard nuclear density profile formula, used in Ref. 
does not apply for A < 40. 

3 He: We followed the parameterization given in Ref. 
|43] ]. which assumed 3 He has a Gaussian density 



profile with core density of 0.93/9,, 
fm. 



and a of 1.11 



2 H: 



It is perhaps not quite meaningful to consider 2 H a 
dense enough nucleus to expect theoretical models 
of medium modifications to be applicable. How- 
ever, several theoretical models have all the right 
ingredients of p medium modifications due to the 
interactions of pions and resonances with the nu- 
cleons and pions in the nucleus and these mod- 
els do reproduce the vacuum spectral shape cor- 
rectly. Furthermore, while the average density of 
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the deuteron is low, recent models utilizing mod- 
ern NN potentials indicate that the deuteron is 
primarily a toroidal structure with a nearly hollow 
core and a maximum density of nearly 2 p nuc at 
r ~ 1 fm. Here, the S and D-wave deuteron wave- 
functions from Ref. 0] using the Argonne AV18 
potential were parameterized and used to calculate 
the AI d = state p(r, 9). 

This density information was used two ways. First, 
it is used to select the appropriate spectral function 
Im Dp(p' B ) for the event. Second, the generated event is 
weighted according to the density at the reaction vertex, 
so that production at the diffuse edge of the nucleus is 
less probable than from the denser core, consistent with 
the requirement of a quasi-free reaction mechanism. 

To determine q, the produced p° is Lorentz- 
transformed to the struck proton plus recoil-fragment 
rest frame, thus specifying the p° propagator value for the 
event, Im D^p^^m^^^q). This value was incorporated 
into the statistical weight for the simulated event. As the 
spectral function varies rapidly with but slowly with 
p' B and q, the propagator was interpolated from the clos- 
est m T7r bin, but no interpolation was performed either 
over p' B or q. The simulated events were then tracked 
through the simulated experimental detectors and ana- 
lyzed in the same manner as the data. 

2. The Post-Leupold-Mosel model 

The PLM model was implemented in a manner similar 
to the RCW model. In this case, 15 in-medium propaga- 
tors were provided |46| in matrix form for nuclear densi- 
ties from 0.1 to 1.5 p nuc . As with the RCW model, longi- 



tudinal, transverse and spin-averaged functions were pro- 
vided, but only the result using the longitudinal spectral 
function is shown here. The propagators were provided 
as a function of p° L energy oj p in the nuclear medium rest 
frame, as well as p° invariant mass, 

A p = A p(p'B, m ^,^p)- 

and was interpolated from the closest m^n bin, but 
no interpolation was done over uj p or p' B . The same 12 C 
density profile and Fermi momentum distributions were 
used as before. 



3. The Saito-Tsushima-Thomas model 

In the STT model, the p° mass is a function of the local 
density at its production point. We use their parameter 
set B for the mean-field value of the cr-meson 

g a a = sip' B + s 2 p'i + S3Pb 
where si = 214.0, s 2 — —44.3 and S3 = 1.9 and g„a is in 
MeV. The p° mass is then given by 

2. , M 8.58 x lO-^MeV^ 1 ), ,. 
m p = m p - -{g a a)[l - (p ff( r)]. 

The complex manner in which the p° mass and width 
are incorporated into the PWA parameterization of Ref. 
|23j| preclude its use here, and so the simpler parameter- 
ization of the p° line-shape of Ref. [2|| is used, instead. 
In all other respects, p° L are randomly generated within 
the nuclear medium via the quasi-free mechanism already 
described. 



